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ABSTRACT 

We determine constraints on the form of axisymmetric toroidal magnetic fields dictated 
by hydrostatic balance in a type II superconducting neutron star with a barotropic 
equation of state. Using Lagrangian perturbation theory, we find the quadrupolar 
distortions due to such fields for various models of neutron stars with type II super- 
conducting and normal regions. We find that the star becomes prolate and can be 
sufficiently distorted to display precession with a period of the order of years. We also 
study the stability of such fields using an energy principle, which allows us to extend 
the stability criteria established by R. J. Tayler for normal conductors to more gen- 
eral media with magnetic free energy that depends on density and magnetic induction, 
such as type II superconductors. We also derive the growth rate and instability condi- 
tions for a specific instability of type II superconductors, first discussed by P. Muzikar, 
C. J. Pethick and P. H. Roberts, using a local analysis based on perturbations around 
a uniform background. 



Key words: 

matter 



Stars: neutron - Magnetic fields - Magnetohydrodynamics - Dense 



1 INTRODUCTION 

Timing residuals varying on timescales of order months to years have been detected in several pulsars, most spectacularly in 
. 5^ PSR B1828-11, where several cycles of nearly periodic variation have been reported (Stairs, Lyne & Shemar 2000; Stairs et al. 
^ : 2003). For PSR B1828-11, the precession period is Pp 500d ^ 4.3 x 10^ s and the spin period is P* « 0.405 s; interpreting 
5^ the long term timing residuals as rigid body precession then implies a stellar distortion e ~ P*/Pp ~ 9-4 x lO^''. Precession 
affects arrival times in two ways (Cordes 1993; Akgiin, Link & Wasserman 2006): (i) Geometrical residuals arise because the 
pulsar beam crosses the plane formed by the angular momentum of the star and the line of sight to the observer at times 
that vary periodically over the precession cycle, (ii) Variations in the angle between the spin and magnetic axes result in a 
periodic variation of the pulsar spindown torque, causing pulse arrival times to vary periodically as well. Precession models 
that combine these two effects describe the data from PSR B1828-11 adequately (Jones & Andersson 2001; Link & Epstein 
2001; Akgiin et al. 2006). 

Problems with these models remain, however. One is the observation by Shaham (1977, 1986) that vortex line pinning 
can prevent long period precession, substituting instead precession with very short periods (of order 10-100 spin periods, 
rather than 10*) that damps out after perhaps 10* cycles, contrary to observations (Sedrakian, Wasserman & Cordes 1999). 
Although Link & Cutler (2002) showed that the precession amplitude in PSR B1828-11 may be large enough to unpin all 
vortex lines in the crystalline stellar crust. Link (2003) argued that the interaction of (magnetized) core superfluid vortex 
lines with the flux tubes in type II superconducting regions would also prevent long period precession. One way out is that 
the core neutrons are not superfluid, an idea that gets some support from comparing theoretical models for cooling neutron 
stars with observations (e.g. Yakovlev & Pethick 2004, and references therein). 

Even if vortex line pinning is not an issue, the required stellar distortion is problematic. Although the rotational distortion 
of a fluid star is substantial, erot ~ Srot/^grav ~ 7 x W~^R^/A4i 4P^ (for uniform density), where i?* = lO^-Re cm and 
M« = IAMiaMq are the radius and mass, and P* is the spin period in seconds, the bulge in a slowly rotating, self-gravitating 
fluid is always axisymmetric about the angular momentum axis, and cannot result in precession. Only the solid crust of a 
neutron star can support distortions that are flxed in the rotating frame of the star, as are needed for precession. However, 
the crust of a neutron star is not very rigid: its shear modulus is only about 0.01 times the crustal pressure. Consequently, 
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e <C Erot if the crustal distortion is "relaxed" at the current rotational frequency of the star (Baym & Pines 1971; Cutler, 
Ushomirsky & Link 2003). For PSR B1828-11, agreement between the observed and calculated precession frequencies would 
require that the crustal deformation be relaxed at a rotation frequency of about 40 Hz, compared with the present frequency 
of about 2.5 Hz (Cutler et al. 2003). 

An alternative explanation for the precession frequency is that it is due to stellar distortions resulting from magnetic 
stresses. The idea that a rotating, magnetic star must precess goes back about fifty years (e.g. Spitzer 1958). If the magnetic 
field and rotational axes are not lined up, then the moment of inertia of the star is the sum of two axisymmetric contributions 
that are misaligned: the rotational distortion, estimated above, and a magnetic distortion of order emag = i?mag/£^grav. In 
such a case, the star will precess about the magnetic axis with a frequency proportional to the magnetic distortion (Mestel & 
Takhar 1972; Mestel et al. 1981; Nittmann & Wood 1981). 

For the typical inferred dipole magnetic fields of neutron stars, the magnetic deformation is far too small, and the 
resulting precession period is far too long: emag ~ 10~^^Si2i?6^iT4 * dipole magnetic field strength B = 10^^Bi2 G. 
However, substantial internal toroidal fields (e.g. B12 ~ 100) could lead to large enough magnetic distortions to account for 
the precession frequency of PSR B1828-11 (e.g. loka 2001; Cutler 2002). 

Larger magnetic deformations could also result from type II superconductivity in the neutron star's core for a given 
magnetic induction strength in the superconductor (e.g. Jones 1975; Easson & Pethick 1977; Cutler 2002; Wasserman 2003). 
In this paper, we shall examine the distortions of a fluid neutron star induced by the enhanced magnetic stresses associated 
with type II superconductivity. Here wc focus on primarily toroidal fields, partly because they are easier to treat, but also 
because they lead to prolate stellar distortions, which the data on PSR B1828-11 seem to favor at least weakly (Wasserman 
2003; Akgun et al. 2006). We will include a weaker poloidal component that can leak into the stellar magnetosphere, as is 
required for the pulsar to be active. Differential rotation within a newborn neutron star most likely amplifies the toroidal 
component of the field (Thompson & Duncan 2001), but stable configurations will require some poloidal field as well (e.g. 
Braithwaitc & Nordlund 2006). We have developed the (more complicated) formalism needed to treat purely poloidal fields 
in a compressible type II superconductor (Akgiin 2007), and will present those calculations elsewhere. 

As a result of ^So pairing via strong interactions, the protons in the interior of a neutron star are expected to form a 
type II superconductor at baxyon number densities between ~ 0.1 — 0.6 fm~^ (e.g. Baym, Pethick & Pines 1969; Baym & 
Pethick 1975; Easson & Pethick 1977; Elgar0y et al. 1996; Jones 2006; Baldo & Schulze 2007). Magnetic fiux penetrates the 
superconducting region in the neutron star in the form of quantized magnetic fiux tubes. Typically, in a neutron star the 
critical field is Hd ~ 10^^ G, and the magnetic induction is i? ~ 10^^ G <C Hd , so the magnetic field is H ki Hd and is 
approximately a function of baryon density (e.g. Easson & Pethick 1977). 

In the neutron star crust, which exists at densities below ~ 2 x 10^*g/cm"^ (Baym, Bethc & Pethick 1971; Lorenz, 
Ravenhall & Pethick 1993), protons are bound in nuclei, and as a result, superconductivity is suppressed. Magnetic stresses 
in a type II superconductor are ~ HB/Aiv « HciB/4n, and consequently will be about Hci/B 10^ times larger than 
those in a normal conductor with the same B, which scale as B^ /8n (Jones 1975; Easson & Pethick 1977). Stresses of this 
magnitude are capable of distorting the neutron star sufficiently to cause precession of the star with a period of the order of a 
year (Cutler 2002; Wasserman 2003). However, we note that hydrostatic equilibrium requires approximate continuity of HB 
throughout the star, so the induction B„ in the normal region is much larger than the induction Bs in the superconducting 
region: B„ oc {HBsY^^ S> Bg. Configurations with large discontinuities in stress are unstable, so it is unrealistic to embed a 
superconducting region with am anomalously large stress inside a star with otherwise much smaller stress. 

The magnetic force in a type II superconductor is inherently different than in a normal conductor. The difference results 
from the fact that the magnetic free energy in a type II superconductor depends both on the magnetic induction, B (or 
equivalently, Wmag = B^ /8n) and on the proton number density, Up. The proton number density is a function of the baryon 
number density, and consequently can be expressed as a function of total mass density, p. A good approximation is to take 
rip <x. p (Easson & Pethick 1977). On the other hand, in a normal conductor the magnetic free energy is a function of magnetic 
induction alone. 

The purpose of this paper is to determine magnetic field configurations in neutron stars with type II superconductors, 
consistent with hydrostatic balance, and assess their stability. We assume that the magnetic deformations are small, which 
enables a perturbative treatment. We neglect rotational deformations, slow fluid motions and associated viscous effects, which 
can be included at a later stage (extending methods laid out by Mestel & Takhar 1972; Mestel et al. 1981; Nittmann & Wood 
1981). With these solutions we can determine the magnetic distortion explicitly (cf. Cutler 2002, who expressed the distortions 
in terms of averages over unspecified field configurations) . 

Assuming (cold nuclear) matter with a barotropic equation of state p{p) imposes significant constraints on the possible 
variation of the magnetic induction B{r,6) in the star. This is because Euler's equation of magnetohydrostatic balance 
requires that the magnetic force per unit mass be a total gradient (a result well known for normal magnetic equilibria; see 
e.g. Prendergast 1956; Monaghan 1965). The fact that H ~ Hci{p) is a function of r alone to lowest order further restricts 
the range of possible B{r,6). With these constraints, we can evaluate the quadrupolar deformation of the star in hydrostatic 
balance (as well as other multipoles, which are uninteresting for precession). In practice, we only calculate these for the 
7 = 2 polytropic equation of state p = np^, where /t is a constant, but the formalism can be applied to any p{p). Moreover, 
although we only present examples for which H = Hci{p), our formalism applies to any magnetic free energy F{p, B), hence 
H = A-KdF/dB = H{p,B). 

Even with the restrictions imposed by hydrostatic balance in a barotropic fluid, and the density dependence of H, many 
possible B{r,9) are permitted, even when we trim the set of solutions by obvious requirements such as regularity. Stability 
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ought to weed out even more possibilities. To examine this question, we use the energy principle that has proved fruitful for 
normal magnetic substances (e.g. Bernstein et al. 1958; Tayler 1973), extended to superconductors in which the magnetic free 
energy (and consequently H) has arbitrary dependencies on p and B. (Roberts 1981 examined this problem for H oc p.) From 
this stability criterion, we show that the most pernicious axisymmetric instability is the interchange instability (just as in 
normal conductors), and we show how the list of candidate field configurations can be winnowed further by requiring immunity 
against it. The interchange instability can be viewed as a magnetic buoyancy mode. Our detailed treatment of perturbations is 
applied specifically to one-component fluids. Buoyancy due to multi-fiuid composition, which arises as a result of the density 
dependence of the number density of charged particles in chemical equilibrium, will introduce new modes (Reisenegger & 
Goldreich 1992), and may change the interchange instability conditions (Ferriere, Zimmer & Blanc 1999, 2001). We postpone 
a complete consideration of these effects to a later paper, but in i]4.3l we argue that stability constraints on the toroidal field 
shape remain the same. 

For non-axisymmetric perturbations, the character of the energy principle is markedly different in the superconducting 
case. From it we find a specific stability criterion for what we will refer to as the Muzikar-Pethick-Roberts (MPR) instability 
first discussed by Muzikar & Pethick (1981) and Roberts (1981), who showed that for sufficiently weak magnetic induction 
B < 10^^ G, the density dependence of H promotes the formation of domains with and without magnetic fiux. From a local 
stability analysis, we show that this instability only acts for m > (non-axisymmetric) modes and only on very small scales 
perpendicular to the field, corresponding to wave numbers ~ 10* /i?*. We estimate the growth time of the instability on 
these scales to be of order 10^ s for typical parameters, i.e. longer than typical Alfven wave crossing times. Although this is a 
distinctive mode associated with type II superconductors, the fact that it only acts on small length scales may cause it to be 
suppressed by small viscous effects. Moreover, since the instability is local it is likely to be present in a rotating star as well. 
Preliminary calculations suggest that while the stability condition is altered by buoyancy, the unstable MPR mode persists 
and has the same growth rate as in a one-component fluid. 

In this treatment, we neglect rotation and internal fluid motions. Our primary goal is to understand the effects of 
the density dependence of the magnetic free energy F{p, B) on equilibrium and stability. This case has been previously 
treated by Roberts (1981), who considered poloidal flelds in a completely type II superconducting star of uniform density 
and magnetic field H (x p. Here we extend these considerations to barotropic equations of state and magnetic fields of the 
form H{p, B) in fiuid stars with type II superconducting shells. We will be concerned with toroidal magnetic fields in this 
paper, deferring the detailed treatment of poloidal fields to future work. We then calculate explicitly the extent of stellar 
deformation due to the magnetic field. Spitzer (1958) and Mestel & Takhar (1972) argued that, to lowest order, the rotational 
and magnetic deformations can be calculated separately. Then, a misalignment in the rotational and magnetic deformations 
leads to precession, as mentioned above. 

In addition to the proton superconductor, there may be a commingled neutron superfiuid in the core of a neutron star. 
If so, the two superfluids are coupled via entrainment. One consequence is that the vortices in the neutron superfluid acquire 
magnetic flux and therefore couple to the magnetic flux tubes in the proton superconductor. This interaction is expected 
to impede precession (Link 2003). The long-term periodicity observed in PSR B1828-11 may require this interaction to be 
of limited scale, perhaps implying that there is no commingling of the two fluids. Moreover, theoretical models for cooling 
neutron stars suggest that there is no compelling observational evidence for core neutron superfluid (Yakovlev & Pethick 
2004). Although gap calculations generally support the existence of a ^So crustal neutron superfluid and a core proton 
superconductor, the theory is less certain about the '^P2 core neutron superfluid. (Elgar0y et al. 1996; Baldo & Schulze 2007). 
Here, we assume that there is no core neutron superfluid overlapping with regions of proton superconductivity. This simplifles 
the problem, as the behavior of a mixed superfluid-superconductor system can be very complex (Glampedakis, Andersson & 
Jones 2007). Moreover, for the reasons given above, this may even be justifled. 

Here, we are primarily concerned with the equilibrium structure of the magnetic field. Although we will also discuss the 
stability from an energy principle point of view, we will not delve into the more comprehensive treatment of modes which 
should also include rotation, internal velocity fields, multi-fluid components, and the elastic crust, as well as dissipation, 
mutual friction and entrainment effects, which would arise in a superfluid-superconductor mixture. In particular, dissipation 
is strongly dependent on whether the neutrons are superfluid or not. Moreover, there will be friction on the magnetic flux tubes 
which is especially important if they coexist with neutron vortices. Stability of rotating stars is known to be affected by normal 
magnetic flelds (Glampedakis & Andersson 2007), and we expect the same to be true in the presence of superconductivity. 
Therefore, our work is only a first step towards a more complete treatment of the neutron star interior, where we highlight 
features arising from the density dependence of the magnetic free energy. 

The outline of this paper is as follows: in ^ we discuss the magnetic stress tensor and force in a type II superconductor. 
In 521 we determine the form of the toroidal magnetic fields in the normal and superconducting regions, consistent with the 
boundary conditions at the stellar surface and internal boundaries. We then proceed with the calculation of the hydrostatic 
equilibrium in the presence of such magnetic fields in various neutron star models with type II and normal regions. We 
calculate the density and gravitational potential perturbations and determine the moments of inertia of the perturbed star. In 
34] we discuss the stability of toroidal fields in the normal and superconducting cases. We show that the interchange instability 
is the worst axisymmetric instability, and derive the MPR instability conditions and relevant time and length scales from a 
local analysis. In 33 we discuss the possibility of adding a small poloidal component to help stabilize the toroidal fields. We 
derive the form of this poloidal field that is consistent with the requirements that the magnetic force be a gradient and that 
the magnetic induction be divergenceless. 
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2 MAGNETIC FORCE IN A TYPE II SUPERCONDUCTOR 

The magnetic stress tensor in a type II superconductor is given as (Easson & Pethick 1977), 



dF dF 



+ ^ . (1) 



The magnetic free energy F{p, B) is a function of mass density, p and magnetic induction, B. In isotropic media the magnetic 
field Hi and induction Bi are parallel, so that Oij — cjji. In general, the relation between the magnetic field and induction is 
given through (Josephson 1966), 

dF 

H = A.-. (2) 

In a normal conducting medium we have H — B, i.e. the magnetic field is independent of density, and the free energy is equal 
to the magnetic energy F — B'^ /Siv. Thus, the stress tensor in this case reduces to, 

B^ BjBj 

On the other hand, the magnetic field in a strongly type II superconducting medium, such as the proton superconductor in 
a neutron star, is H ~ Hci 2> B, and depends most sensitively on the proton number density Up and the superconducting 
energy gap A (Tinkham 1975; Easson & Pethick 1977), which are functions of baryon density p (Elgar0y et al. 1996; Baldo 
& Schulze 2007); therefore, H ~ H{p) and F ~ HB/A-k. In this case, the magnetic stress tensor reduces to, 

<x.=-p^5.. + ^. (4) 

The stress tensor used by Roberts (1981) is of this form, with H oc p. 
In general, the gradient of the free energy is given as, 

dF dF 
From equation ((2} it follows that. 

Making use of these relations as well as the fact that V ■ B = 0, the magnetic force density can be calculated from equation 
Q as, 

r yj ^ dF dF ^ B.VjH, [(V xH)x B], ^ dF 
- V,-. = -pV.— - BV.— + ^— = pV.- . (7) 

This is the form of the force in a type II superconductor. (In fact, it is true in any magnetic medium where the free energy 
is a function of density and magnetic induction.) This is inherently different from the force in a normal conducting medium, 
which can be retrieved by setting H — B and F = B'^ 
In hydrostatic balance, 

Vp + pV(^ = /^,g , (8) 

where p is pressure, p is mass density, (j> is gravitational potential, and f-^^^ is the magnetic force density (equation [T]). In 
barotropic equations of state, pressure is a function of density and we can define dh{p) = p~^dp{p); then, 

pVih + ^) = f^^^. (9) 

This equation requires the magnetic force per unit mass to be a gradient of a potential, i.e. f^^^^ = —p^il). We will express 
the magnetic potential as the sum of two terms, 

= + , (10) 

where, we define, 

{V xH)xB JxB . . dF 

^ = = -pVi/)i and tpu = ^ ■ (11) 

47r c dp 

J is the current density, ipi is the magnetic potential for a normal conductor, and ipu is present only for a type II superconduc- 
tor. The second term in the magnetic force (equation [T]) is already a gradient. On the other hand, note that the requirement 
for the first term to be a gradient can be expressed alternatively as, 

Vx(i^).„. 

This equation needs to be satisfied for both the normal and type II superconducting cases, and imposes a severe restriction 
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Figure 1. A star with a normal core, superconducting sliell, and a surrounding normal layer. The radius of the core is ri and the outer 
radius of the superconducting shell is r2. 

on the form of the magnetic fields, which are also required to satisfy V • B = 0. The normal conducting case is discussed, for 
example, in Prendergast (1956) and Monaghan (1965). For the strongly type 11 case and H (x p, Roberts (1981) found poloidal 
field configurations for uniformly dense stars, and Akgiin (2007) found poloidal field configurations for 7 = 2 polytropes. 



3 TOROIDAL FIELDS 

The current density for a toroidal field H = H{r,9)4> is, 

— = V X = V(//rsine») X — ^ . (13) 
c r sin 



Taking the induction to be B = B{r,6)(j), we get, 
JxB_{VxH)xB_ BViHrsmd) 



(14) 



pc Amp Aixpr sin 6 

This is clearly a total gradient, as required by equation (|12p . for magnetic inductions of the form, 

B{r,6) ^ 4npr sm6f{Hrsm9) , (15) 

where / is an arbitrary function of ^ = HrsinO. The factor of 4-k is included so that defining a new function through 
/(C) = g'iO gives, using the definitions in equation (fTTjl . 

l2i^^-Vg{0 i.e. Mr,0)=g{(;). (16) 
pc 

This is valid for any H{r,8). However, for a strongly type II superconductor H ^ H(r), and we have (equation llip. 

B dH d\nH ^ ,,„, 

^11 = ■7--r = -71 — (9 (0 ■ 17 

47r dp dinp 

For a normal conductor H = B, and equation (|15|l implies that the magnetic induction is now given through the form, 

Bir,9) = t<£I^^, (18) 
r sm 6^ 

where h is an arbitrary function of ^ = pr^ sin^ 6. It then follows that, 

^^-V^^-^ i.e. ^'(C)^^^. (19) 
pc 87rf 47r^ 

Note that, for a uniform density, the magnetic induction is a function of the cylindrical radius, ru — rsinO. 



3.1 Star with a Superconducting Shell 

Consider the case of a strongly type II superconducting region confined to a spherical shell between radii ri and r2 (where 
r2 > ri). Let the magnetic field be Be inside the normal core, H inside the superconducting shell (with a corresponding 



6 T. Akgiin and I. Wasserman 



magnetic induction Bs), and B„ inside the normal outer layer (as depicted in fig.[T|. Since the fields have no radial components 
in this case, they need not be continuous across the boundaries, and there will be surface currents. 

In fact, it turns out that in the toroidal case it is not possible to have a continuous magnetic field across the boundaries, 
if H — H{p) in the superconducting region. Consider one of the boundaries of the superconducting shell, located at r = rj,. 
For the present discussion, it is immaterial whether the normal region lies on the inside or the outside of the boundary. In 
the absence of surface currents, the boundary condition that follows from Maxwell's equations requires the continuity of the 
tangential magnetic field, 

f X H ^ f X B„ . (20) 

Since H is & function of radius in a strongly type II superconductor, for this equation to be satisfied everywhere on the surface 
of a spherical boundary, the magnetic field _B„ inside the normal region (given by equation [18} would have to be a function 
of only radius at the boundary as well. This implies that we must choose a function h{^) oc £,^^^, so that Bn{r,6) oc p^^'^(r). 
However, in this case, the magnetic potential becomes ipniO oc In^ (equation I19|) . which diverges whenever ^ = pr^ sin^ ^ is 
zero. In other words, it diverges at the center of the star {r —> 0), at the surface (p — > 0), and along the symmetry axis {8 —> 0). 
We also note that when the magnetic induction Bs inside the superconducting region (given by equation I15|l is chosen so that 
it is angle independent (i.e. f{() oc the corresponding potential is also logarithmic, 'i/)i(C) oc ln(". 

We therefore conclude that continuous toroidal fields, or more generally, angle-independent magnetic inductions, are 
inconsistent under the assumption that H = H{p) holds up to the boundaries of the superconducting region. In a more 
realistic treatment, H{p,B) should be allowed to decrease smoothly to about Bs near the boundaries, which would remove 
the need for surface currents. 



3.2 Boundary Conditions 

Hydrostatic equilibrium for a fluid with a barotropic equation of state, in the absence of magnetic fields, is spherically 
symmetric and is given by (from equation O, 

V{h + <p)^0. (21) 

When a magnetic force that is small in comparison to pressure and gravity is applied, the equilibrium quantities are changed 
by small amounts 5p, Sp, 5h and 5<j), where S denotes Eulerian changes. Writing the magnetic force in terms of the magnetic 
potential, f^^^ = —p^ip, the equation for the perturbations around the background equilibrium can be written as, 

V{Sh + 5(j> + ilj) =0 . (22) 

From here it follows that, 

5h=^Sp^^o-S(l>-i^ . (23) 
dp 

58o is Bernoulli's constant and is the same for the entire star. This can be understood by treating the entire star as a single 
fluid region, with a magnetic potential that varies continuously throughout the interior, but that has steep changes in some 
small intervals corresponding to the boundaries. 

While the background quantities p, p and (j) are continuous throughout the star, their perturbations are not. Only Sij) 
and its gradient are required to be coiitiimous, since there cannot be delta functions in mass. This implies that there will be 
a density perturbation jump at a boundary, given by (from equation [23]) , 

^((Sps - Spn) = -tps + tpn ■ (24) 

dp 

Here the subscripts s and n refer to the superconducting and normal regions, respectively. 

There must be substantial surface currents at the boundaries of the superconducting shell, and therefore, the magnetic 
fleld is discontinuous across them. Otherwise, as discussed before, the magnetic potentials become singular. From the continuity 
of stress, it follows that, 

(25) 

Eij is the total stress tensor and Uj is the normal unit vector of the boundary, which in this case is simply the radial unit 
vector f. Thus, we require the rr, rO and rcj) components of the stress tensor to be continuous. The last two vanish identically 
for fluids with toroidal flelds. 
The total stress is, 

Sij = —SpSij + (Jij , (26) 
and from equation (|25p . we have, 

— Sps + (Jrr.s ~ —Spn + Orr.n ■ (27) 

Using the fact that for a polytrope p — Kp'^ , we have dh/dp = jp/p'^ and 5p — {'yp/p)Sp, we can combine this result with 
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equation (|24|) to get, 
TP 

{Sps — 5pn) = —p{lps — Ipn) = CTrr.s " CTrr-.n ■ (28) 

P 

The components of the stress tensor inside the normal and superconducting regions are given by (equations |3] and |4| , 

O-rr.n = — — aud (T^r.s = ^ = ^ • (29) 

Using t/js = i/ji + t/jii (equation [T0)>. we thus obtain, 

= + (30) 

This equation needs to be satisfied by the magnetic fields at the boundary. Note that since cx HBs/p and ipn oc B^/p, 
this equation implies that B„ oc (HBg)^^'^. If we take H ^ Bg to hold at the boundaries of the superconductor as well as its 
interior, then the boundary condition clearly requires B„ 3> Bs- Taking a more general H{p,B), varying continuously from 
-ffci (p) to Bs through a thin boundary layer, would result in a smooth but similar growth in the magnetic induction between 
the strongly type II and normal regions. (Surface currents would be smoothed out over this boundary layer.) For entirely 
normal conductors, the corresponding boundary condition simply implies the continuity of magnetic fields. 

In a more sophisticated treatment of the transitions from superconducting to normal and/or fluid to crust, two dimen- 
sionless ratios characterize the superconducting state. One is, 

-7«¥^- (31) 

where A is the London penetration depth, ^ is the coherence length in the proton superconductor, Up = 10''^np.37 cm~^ is the 
proton number density, and A is the proton superconducting gap. The other is, 

j^68B;,'^'{np,sr)'^\ (32) 

where a is the spacing between flux tubes (Tinkham 1975). In a type II superconductor, k > l/\/2- 

At the crust-core boundary, Up falls dramatically, and a/A drops, which means that interactions between flux tubes 
become important. As a result, our approximation that H ~ Hc\{p) must fail, and must be replaced by a more general (and 
complicated) function of both p and B. 

At the inner boundary of the superconducting layer, A ultimately disappears, and k falls below l/\/2- In this regime, we 
expect a boundary layer of a type I superconductor to form. In fact, it is also possible for such a layer to form at the crust-core 
boundary, since the gap depends exponentially on the density of states near the proton Fermi surface, which falls with proton 
density. Thus, at both boundaries, we expect the magnetic fleld to decrease rapidly from H ~ 10^^ G to Bn ~ (HBs)^^^ ■ 



3.3 Derivation of the Magnetic Fields 

We will assume a simple power law relation between the magnetic fleld in the superconducting region and mass density, 

(33) 

where He and pc stand for the central values of the corresponding quantities. When the superconducting region is conflned 
to a shell, we can take He to be the extrapolated fleld strength at the center. In reality, in a strongly type II superconductor, 
H depends on the superconducting energy gap A, in addition to the proton number density Up (Tinkham 1975; Easson & 
Pethick 1977). Both Up and A are functions of baryon density p (Elgar0y et al. 1996; Baldo & Schulze 2007). A vanishes at 
sufficiently high densities, and protons become normal. At low densities, superconductivity is suppressed since protons are 
bound in the nuclei in the neutron star crust. In both cases, the transition from superconducting to normal state may be 
sharp and we take the form given by equation (|33p in superconducting regions. 

In this case, equations (|16p and (|17|) imply ipi = g{(^) and ipu = aCd'iO^ where = HrsmO. Consider a power law 
function of the form g[Q — N(", where A*' is a constant; then tpi — N(" and iIju = naN(", so that the total magnetic 
potential becomes, 

i's=ipi + i^u = {na + 1)NC . (34) 

We exclude n = since that corresponds to zero magnetic induction and force. On the other hand, for n < the magnetic 
potential diverges when either r — > or S — » 0. Moreover, the magnetic force diverges in the same limits in the interval 
< n < 1 . Therefore, the only nonsingular choices are n ^ 1. The magnetic induction inside the superconductor is (equation 

\m, 

/ „ \ CT(n-l) + l / \ " 

Bs{r,9) = Boiy] i—j sm"e where = 47rnA'pci/ ""Vo" . (35) 
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The constant Vo will be defined later. The corresponding magnetic potential can be written as, 

(ncr + l)H,Bo 



71(7 

p \ I r 



where 



Pc J \ro J 'innpc 
Inside the normal region we have, from equations (|18|) and (|19|l . defining f — pr^ sin'^ 6, 



Bn(r, 



and ip'niO 



r sm y . 

We will assume a power law for the arbitrary function, /i(^) = M^™, where M is a constant. Then, 



and ijjn ~ 



8np Stt " 47r(2m - 1) ' 

The boundary condition (equation I30|l gives, after some rearrangement. 



ivc" 



(36) 



(37) 



(38) 



(39) 



87r(2m - 1) 

In order to satisfy this equation for all values of 6 at the boundary (which we will assume to be located at some radius r — rt 
we must have, 



n = 4m — 2 whence M 



A-KnNH"{n) 



P'^/^in) 

Then the magnetic field in the normal region is 

/ s (n + 2)/4 / „ \ n/2 



1/2 



where Bo = M p, 



(n + 2)/4 n/2 



(40) 



(41) 



Note that Bs and B„ must have different angular dependencies in order for the potentials tps and il^n to be consistent. 
Moreover, 



Bo ~ (HcBo 



1/2 



pin) 



-1 n(2(T-l)/4 



(42) 



so that the magnetic fields in the normal regions are moderately strong. The magnetic potential in the normal region is, 

(n + 2)Bl 



P \ I r 



\pc/ \ro 
Thus, it follows that *o oc *o. 



where '^o ~ 



STvnpc 



(43) 



n + 2 Bi 



n + 2 



pin) 



-1 n(2cT-l)/2 



*o 2{na + 1) HcBo ^ 2(na + 1) 
As in the superconducting case, we need to have n 1 in order to avoid any divergences in the potentials or forces. 



(44) 



3.4 The n = 1 Case 

In a later section, we will show that toroidal fields by themselves are unstable, and that the n — 1 case is the closest to being 
stable. We will be concerned particularly with cases where H oc p, i.e. a — 1. This corresponds to taking the proton number 
density to be proportional to the baryon density, Up oc p, and neglecting logarithmic dependencies in H , which is a good first 
order approximation (Easson & Pethick 1977; Muzikar & Pethick 1981). The magnetic potentials in the superconducting and 
normal regions become, from equations (|36p and (I43p . 

(45) 



^ *° (p~) (^) ^^"^^ and ^" = *° (^) {y^^^^^ 
where, from equation (|44p . we have. 



HcBo 

2-Kpc 



, *o 3 
and -;— = T 

4 



Pi^b) 



Pc 



1/2 



The angular part of the potentials can be expanded in Legendre polynomials, 

oo 

sin 6* = ^QiPpXcose) . 



(46) 



(47) 



Only even I remain in the series and the coefficients are. 



2/ + 1 

Qi.= ^ sin6»Pf(cos6')d(cos6l) = 



(21 + l)7r2 



2{l + 2)(1 - £)r2(^/2 + l)r2(l/2 - £/2) 



(48) 
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Figure 2. Magnetic potential profile for a three component star witfi a normal core, type II superconducting shell, and surrounding 
normal layer. The potential is shown for the n = a = I case for the magnetic field (equation I45I I. and a 7 = 2 polytropic equation of 
state. The superconducting shell lies between xi = O.Gtt and X2 = O.Qtt, and is shown shaded. The potential is shown along the equator 
of the star, i.e. sind = 1, in units of >I'o defined in equation 1 146 I I. The profiles for the potentials within each region are shown extended 
over the whole star for comparison. 



In particular, Oo = 
Qe+2 ^ i2£ + 5)i£- 



r/4. Subsequent terms in the expansion have the ratio, 

i){e-i) 



(2^ + l)(^ + 4)(£ + 2) 



(49) 



Clearly, Oi+2/Oe ^ 1 as £ —> <x. The result can also be expressed in terms of the spherical harmonics which are related to 
the Legendre polynomials through. 



2^ + 1 



4-K 

Then, for even I 



Pi (cos t 



sine = ^ QeYeie) where Qe 



47r 
2ITT 



(50) 



(51) 



We will consider a 7 = 2 polytrope for which the equation of state is p = np^, where k is a constant. In this case, the 
background density is of the form p = pc sin a;/a:, in terms of the dimensionless variable x = r/vo, where Vo — \J k/2tiG. 
The stellar radius is J?* — nro, and the stellar mass is M* = ttMo, where Mo = Anpcrl. The central density is given by 
pc = nMi^/ARl. For a neutron star with Af* I.4M0 and iJ* ^ lO'' cm, we have pc « 2.2 x lO^'^ g/cm^. 

As noted before, superconductivity exists only within a certain range of densities, or equivalently, a range of radii, which 
we will denote by a;i < a; < X2- In particular, it is suppressed in the crust where the protons become bound in nuclei. The 
crust exists at densities below p « 2 x 10^'*g/cm^ (Baym et al. 1971; Lorenz et al. 1993), corresponding to an outer radius 
of X2 ^ 0.97r. On the other hand, the proton pairing gap vanishes at higher densities. This cutoff for superconductivity is not 
as well-established and estimates ranee from p « 5 X lO" g/cm^ to 10^^ g/cm^ (Elgar0y et al. 1996; Baldo & Schulze 2007). 
Thus, the inner boundary of the superconducting shell ranges from xi ~ O.Stt to O.Bvr, respectively. 

The magnetic potential for the n = 1 case in a three component star consisting of a type II superconducting shell 
surrounded by normal regions (as depicted in fig. [l| is shown in fig. [2] Note that the potential within the superconducting 
shell (which is taken to be in the interval 0.6 < x/n < 0.9) is larger than those in the normal regions. 



3.5 Calculation of the Gravitational Potential Perturbation 

The gravitational potential perturbations are given by the perturbed Poisson's equation, 
V^(5<^ = AttGSp . 



(52) 
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For a 7 = 2 polytrope, we have dh/dp = p'{p)/p = 2k, and equation (|23|) becomes 2k(5/9 = 23o — 54' — lii. Expanding the 
perturbations in spherical harmonics as S(j){x,0) — (pti{x)Yi{9) and so on, Poisson's equation gives. 



d 



^2 

x^ dx \ dx 



(f)i = ^oSio — tpe ■ (53) 



The complete solution of this equation is the sum of a homogeneous solution and a particular solution. The homogeneous 
solution is given in terms of the spherical Bessel functions, 4>h{x) — Aeji{x) + Biyi{x), and the particular solution can be found 
by the method of variation of parameters, <j>p{x) = Ai{x)ji(x) + Bi{x)yi{x). Thus, the gravitational potential perturbations 
in the three regions (core, superconducting shell and outer normal layer, as depicted in fig. [TJ are, 

(l>c,e{x) = [Ai + Ae{x)^ je{x) + ^Be + Be{x)^ ye{x) + ^oSeo 

<l>sA^) = [Ct + Ciix)'^ Jt{x) + [Oe + De{x)j ye{x) + ^^Sm (54) 
(l>„,e{x) = ^Ee + Ee{x)'^ ji{x) + \Ft + Fi{x)^ yi{x) + ^oSm 
where Ai through Fi are constants, and we define, 

Ai{x) = - / t^iP,At)ye{t)dt and Be{x) = - / t^ij,,e{t)je{t)dt . (55) 



J X Jo 

Here xpc,e refers to the f-th component of the spherical harmonic expansion of the potential t/jc. The remaining coefficients 
are defined in an analogous fashion. Note that the integration boundaries can be arbitrarily adjusted, which amounts to a 
redefinition of the constants Ae through Fe above. The particular choice made here makes sure there are no singularities, but 
is otherwise immaterial. 

Since there can be no gravitational forces in the center, the gradient of the gravitational potential must vanish there. 
This implies that as a; —> we must have (l>e ^ constant for ^ = 0, and — > and (fi'i^O for ^ 7^ 0. As a; — > 0, the limiting 
values of the spherical Bessel functions are ji oc x^ and ye oc x~^~^. It therefore follows that B^ = for all values of £. The 
remaining five coefficients At, Ci, Di, Ei and Ft, and Bernoulli's constant *8o are to be determined from the continuity of the 
potentials and their derivatives across the shell boundaries, which we will take to be located at xi and X2, such that x\ < X2, 

4>c,e{xi) = (ps,i{xi) and (p'i{xi) = (^'3 f(xi) 

(56) 

(psji^i) = (I>n,i{x2) and <?f>l,f (2:2) = (I>'„,i{x2) 
and from the boundary conditions at the stellar surface, which is located at x = tt, 
TT^'n.eM + l)4>n,t(n) =0 for ^ / 

(57) 

0n,£(7!-) = 0n,«(7r) =0 for ^ = 

The surface boundary conditions follow from the multipole expansion of the gravitational potential, which implies that 
(f>i oc x~^~^, and the conservation of mass, which additionally implies (pi — for £ = 0. 

Making use of various relations between spherical Bessel functions^ the continuity conditions at the shell boundaries 
(equation I56|) yield, 

Ae + Ae{xi) ^ Ci + Ce{xi) and Be(xi) ^ De + Diixi) 

(58) 

Ce + Ce{x2) ^ Ee + Ee{x2) and De + De(x2) = Fe + Fe(x2) 

and the surface boundary conditions (equation I57|l give, since Eilji) = 0, 

Etji^ii-jx) + \Fi + F,(7r)] yi-i{TT) = for l^Q 

Ee Fe + FeM ^^^^ 

Tv^yi(Tv) 7r^ji(7r) 

Special cases can be considered. For instance, for xi = and X2 ~ it we retrieve the completely superconducting star. In 
this case Be{xi) — De(xi) = so that De = 0. Since Ce{x2) = as well, the surface boundary conditions reduce to, 

Ceje~i(n) + De{n)ye^^{n) =0 for £ / 

^-^-4^ for £.0 ^''^ 

On the other hand, letting xi — ^ while keeping a;2 < vr we retrieve the case of a superconducting core surrounded by a normal 

In particular, letting fe denote either je or ye, we have je{^)yl {x) — je (^)yei^) = , = xfe—i{x) — (I + l)fe{x) and 

(21 + l)fe'(x) = efe-i(x) - {£ + l)fe+iix). 
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Figure 3. Gravitational potential perturbation for a fluid star with toroidal fields, expanded in spherical harmonics for the n = 1 case 
(equation I45| l. The potentials are shown for four sample models: type II superconducting shell between x\ = O.Stt and X2 = O.Qtt (case 
A) and between xi = O.Gir and X2 = O.Qtt (case B), completely superconducting star (xi = and X2 = tt), and completely normal star 
(3:1 = X2 = O.Qtt). The figure on the left shows the first two harmonics 0^ (for £ = and £ = 2) scaled by the maximum value of the 
potential, flimax. The figure on the right shows </<max for the first few £, in units of defined in equation II46I I. The points for different 
values of £ (shown with circles) are connected by a cubic spline curve. The amplitude of (pi decreases sharply with £. 
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Figure 4. Gravitational potential perturbation for a fluid star as a function of £. The potential is shown for the n = 1 case of a three 
component star with a superconducting shell between xi = O.Gtt and X2 = O.Qtt. The same scaling is used as in fig. [3] and only £ > are 
shown. 

region. When xi = X2 the star is completely normal conducting. All such cases are equivalent, up to a scaling determined by 
the magnitude of the magnetic potential (which is given through equation [44]) . Sample models are shown in figs. [3] and |4] for 
the n — 1 case discussed before (equation I45p . 

3.6 Density Perturbation 

The density perturbation within each region can be calculated through equation (|23[) . which for a 7 = 2 polytrope becomes, 

2k5p = ®o - (5<?!> - V • (61) 

Sample plots of density perturbations for the n — 1 case are shown in fig. [5] The density jump at a boundary is then given 
through, 



2k Ap = 2K{5pin - Spout) = V'out - V'in 



(62) 
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Figure 5. Density perturbations for a fluid star witli toroidal fields, expanded in spherical liarmonics. Plots are shown for the same four 
sample cases considered in fig. |3] The figure on the left shows the first two harmonics pi (for £ = and £ = 2) scaled by the surface value 
of the density perturbation (tt). The shaded regions indicate the position of the superconducting shell. The figure on the right shows 
2Kpe{iT) for the first few £, in units of '^o defined in equation ((46). 



In particular, consider the density jump when going from a normal region into a superconducting region at a boundary r = rt- 
Using equations (|36|l and (|43|) . we get, 

2kAp = 2k((5p„ - Sps) = ips ~ ipn ^ oP"^ , ^l i's{rb,d) where n ^ 1 . (63) 

Zyna + ij 

Note that Ap ^ for o > 1/2. In other words, the density perturbation decreases when going from a normal region into a 
superconducting region, and vice versa. Also note that the jump goes to zero at the poles, i.e. Ap ^ as ^ ^ 0, since the 
magnetic potentials vanish there. 

The relation between the Eulerian density perturbation and the Lagrangian displacement is given through, 

5p = -V- (pi) = -pV ■ { - p'^r . (64) 

Normally, the term V ■ { inside the fluid is undetermined. However, at the surface p = 0, so that we can calculate the radial 
displacement, which determines the shape of the perturbed stellar surface, 

= -Sp/p' . (65) 

For a 7 = 2 polytrope we have p = pcsinx/x, so that at the surface p'(7r) = —pc/n and = ivSp/pc. The £ — term in the 
spherical harmonic expansion of defines a spherically symmetric expansion (or compression) of the star, while higher order 
£ determine the deformation of the surface as a function of the polar angle, 8. 



3.7 Quadrupolar Distortion 

The moment of inertia of the unperturbed star is given by, 

I^j = p{r'^5^J -r^rj)(fr . (66) 
Jv 

Since the star is initially spherically symmetric we have Ixx = lyy = Izz- For a 7 = 2 polytrope the density profile is given 
through p = pc smx/x, so that the moment of inertia becomes, 

r - r f 2 /, ■ 2 a 2 \ ,3 8(71^ - Q)pcRl .„„s 
lo = Ixx = / pr (1-sm fc'cos ifi) d r = — . (67) 

Here ii* is the stellar radius, which corresponds to a; = Ri,/ro ~ it. 

The application of the magnetic perturbation renders the star axisymmetric (Ji = I2 7^ la)- In this case the moments of 
inertia become h = lo + &I1 around an axis that lies in the equatorial plane, and Is = Jo + SIz around the axis of symmetry 
which passes through the poles. We will define the star to be oblate when 5/3 > 5Ii and prolate when Sis < SIi. In other 
words, when more of the mass is distributed towards the equator the star is oblate, and when more of the mass is closer to 
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Figure 6. (f>2{R*) as a function of tiie width w = {x2 — Xi)/x2 of the superconducting shell in a three component star. The outer radius 
of the shell is fixed at X2 = O.Qtt. The type II shell models (cases A and B) listed in table[T]are retrieved by setting xi = O.Sir (w = 1/9) 
and xi = O.Gtt {w = 1/3), respectively. When xi = X2 {w = 0) the star becomes normal. 



Table 1. Values of (f>2{R*) for the cases considered in fig. [3] The negative signs signify the fact that the models considered here are 
prolate, i.e. SIi > SI3. 



Case 02(R*)/*o 



type II shell (A) 


-1.67 


type II shell (B) 


-2.18 


superconducting 


-2.33 


normal 


-1.18 



the poles the star is prolate. The diflterence between the moments of inertia is related to the gravitational quadrupole moment, 
which in turn is related to the £ = 2 harmonic of the gravitational potential at the stellar surface, 

Q20 = l^pr'Y2{9)d'r = ~^{5h - Sh) = -^-^^^ ■ (68) 
Thus, 

[^ G{5h-5h) 

MR*) = y ^ ^3 • (69) 

Therefore, the sign of 02 at the surface determines whether the star is prolate or oblate. Note that for all the cases shown in 
fig- [31 4'2{R*) is negative and consequently the star is prolate. The precession frequency of an axisymmetric star is ~ ef2*, 
where f2* is the angular velocity and e is a dimensionless constant defined through. 



h " Jo 16(7r2 - 6)Gpci?2 



(70) 



For the n = 1 case, the gravitational potential perturbations are measured in units of = HcBo/2-Kpc (equation [46]). The 
central density for a 7 = 2 polytrope is pc — -kM^ / A:R^. Thus, we can rewrite the above equation as, 

, = 945 X 10- ("MM^ (Jl^\ " " (71) 

Sample values of (j)2{R*) are listed in table[l] and (j}2{R*) as a function of superconducting shell width in a three component 
star is plotted in fig.|6] Note that the values of e for the various models are very similar. This should not be surprising, as the 
magnetic fields in all cases are of similar magnitude. 

In particular, the normal case considered here (in figs. [3] and [5l and in table [1} is for a magnetic field of strength 
Bo ~ {HcBoY^'^ [p{^2) / pc]^^'^ ~ 1.8 X lO^'' G (equation I42|l . This is simply the limiting value of the normal field as the 
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superconducting shell vanishes, xi X2. In the normal case, the magnetic potential is given in units of '^o = SB^/Sirpc 
f equation 143 [). which can be evaluated for different choices of _Bo- 



4 STABILITY OF MAGNETIC FIELDS 

In this section, we will discuss the stability of toroidal fields in neutron stars. We will follow the energy principle considerations 
outlined in Bernstein et al. (1958) and Tayler (1973). An extensive review is also given in Freidberg (1982). The formalism 
that is developed in this section is valid for any H{p, B) and is applicable to both normal and superconducting neutron stars. 
For the purpose of this section, we will treat the entire star as either normal or superconducting, and therefore will not worry 
about internal boundaries. 

We also ignore rotation, and thus do not need to pay attention to "trivial" displacements discussed by Friedman & Schutz 
(1978). In magnetic stars, trivial modes are defined by the requirements that 5p = Q and SB — 0. Since we will express the 
energy of the perturbations in terms of Sp and SB, trivial displacements will have no effect on it (see equation B60 in Friedman 
& Schutz 1978 and footnote 3 in Glampedakis & Andersson 2007). However, in a rotating star, trivial displacements will have 
to be taken into consideration. 

Glampedakis & Andersson (2007) emphasize the importance of the magnetic field for rotating stars by showing that 
sufficiently strong fields can stabilize inertial modes that would otherwise be unstable. The same will be true for type II 
superconducting stars. We will not treat rotation-induced instabilities here. Instead, we emphasize the effects of the magnetic 
free energy F{p, B) in a type II superconductor. Energy conditions presume zero dissipation. Moreover, we consider a single 
fluid, which in reality consists of at least three fluids: neutrons, protons and electrons. There will be additional buoyant modes 
which may or may not alter the stability conditions we derive. 

Assuming small oscillatory perturbations about equilibrium, we have, from equation (|8]), 

- P§ = pu\ = s{Vp + pV4> - f^^J = -:f{^) . (72) 

The force operator is self-adjoint, which implies that the eigenvalues uj'^ are real. One condition for stability is that all 
frequencies u be real, so that there are no growing modes. Alternatively, the variation in the total potential energy due to the 
perturbations should always be positive, 

SW^-^li-J'ii)dV>0. (73) 

To lowest order, the integration is carried over the equilibrium volume. The Lagrangian and Eulerian pressure perturbations 
are given by Ap — {'yp/p)Ap — — 7pV ■ ^ and Sp = Ap — ^ ■ Vp = — 7pV ■ ^ — ^ • Vp. Here 7 is for the perturbations, and 
in general may differ from the background polytropic index. The difference gives rise to buoyancy terms, which will not be 
considered in this paper, however we will comment on their effects on stability briefly. 
Integrating by parts, we get, 

SW = SWp + (JW'mag 

SWp = If [7P(V ■ + (I ■ Vp)(V ■ I) - (I ■ V,^)(V ■ p|) + Pi ■ dV 

1 / r 1 (^4) 

--^ dS-l [7PV-C + C- VpJ 



We will refer to the two parts in the energy as the hydrostatic part SWp, which includes the contributions from pressure and 
gravity, and the magnetic part SlVmag- In equilibrium, the pressure and density are related through a polytropic equation of 
state and consequently they both go to zero at the surface. Therefore, the surface integral vanishes. 

We now turn our attention to the calculation of the magnetic energy variation. Faraday's law gives the variation in the 
magnetic fleld in a perfect conductor as, 

= V X ({ X B) . (75) 

We next discuss the normal and superconducting cases separately. 



4.1 Normal Conducting Star 

In a normal conducting medium, the force is given as, 
^ JxB ^ {VxB)xB 

•/mag ^ 4^ • I'"; 
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The perturbed force becomes, 



SJ X B J xSB 

•^JmaE = ^ 



where 



c 



V X SB 



c c c An 

Integrating the first term in SWma^, given through equation (|74|l . by parts and rearranging, we thus have. 



mag 2 



\5B\- 



J 5Bxi 



dV 



8-K 



dS 



■ SB) - B{i ■ SB] 



(77) 



(78) 



The first surface integral vanishes when dS ■ B = 0, i.e. when the magnetic field is perpendicular to the surface, as is the 
case for a toroidal field. On the other hand, the second surface integral vanishes when the field vanishes at the surface. This 
form of the energy variation is the same as that given by Bernstein et al. (1958) for dS ■ B = 0. The surface integrals may be 
relevant, for instance, in the case of poloidal fields. However, we will not need to worry about these as we will be considering 
toroidal fields that vanish at the surface. 



4.2 Type II Superconducting Star 

The magnetic force for a type II superconductor is given by equation ((7}, 

/mag = - PV^II , (79) 

where ^pu = dF/dp, from equation (|lip . The current density is now given through AnJ /c — ^ x H . The magnetic free energy 
_F is a function of p and B and is related to the magnetic field through equation H = AirdF/dB. The perturbation of the 
force gives, 

SJ X B J X SB ^ „ , „„, 

= + SpVi^U - pVSxPu . (80) 

Consider the energy due to the first term of the magnetic force. Following the same procedure as in the derivation of equation 
ifTS]! . we get. 



^Ji-SJxBdV = J i X B ■ {V X SH)dV 



= ^jdS- [b({ ■ SH) - i{B ■ 5if)] - ^ j SH -SBdV 



(81) 



When B vanishes on the surface we can drop the surface integral. On the other hand, note that we can rewrite the last two 
terms in the magnetic energy variation as, 

{Spi ■ V^ii + pi ■ V(5Vii) dV = j {SpS, ■ VV'ii + SpStPu) dV = j (5pAi/>ii dV . (82) 

Here, we have made use of the relation A = 5 + ^ • V, between Lagrangian and Eulerian perturbations. Thus, the magnetic 
energy variation for a type II superconductor becomes, from equation (|74[) . 



U\SJLlE^£llE2il+5pA^u] dV+^idS- kiB-SH)-B{i-SH) 
2 J 47r c 87r J L 



(83) 



This is to be contrasted with the magnetic energy for the normal case given by equation (|78p . In particular, the first two 
terms in the volume integrals are of the same form, with a B in the normal case replaced by an H in the superconducting 
case. The same is true for the surface integral terms. However, in the superconducting case there is also an additional term 
that arises from the potential tpu, that has no analogue in the normal case. 

In the strongly type II superconducting case the magnetic field is a function of density only, H — H{p). On the other 
hand, in the normal case we have H — B. In general, H, -i/jn and F will all be functions of p and B. Using the definition of 
the potential tjju from equation pip , we get, 

d^F d'^F 



We will assume that the form of SB given through equation (|75p is still valid for the superconducting case. Also note the 
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following relations which will be of use, 
SB 5BB 



B B 
5B = B-SB 

SH = 5HB + H5B 

^ dH - OH - ^ 
SH^ —Sp+—SB 
op oB 

Note that B 1. SB, which also follows from S{B ■ B) = 0. Using the above relations we have, 



SH SB = SHSB + ^ 



SB SB [SBY 



Using equation ([2} which relates H and F, the perturbation in the magnetic field can be written as, 
\dpdB ''^ dB^' 



SH = 4n( -^^Sp + 



This allows us to express the energy in terms of derivatives of F. 

For a strongly type II superconductor H oc p, equation (|83p reduces to (Roberts 1981; Akgiin 2007), 



SW^ 



^ J ^SH ■ SB ~ SB ■ i X {V X H) ~ {H ■ SB){'V ■^)+SH -i^- VB) - SB-{^- Vi?)] dV 
^^{SH B + H SB) - B{€, ■ 5if)j 



+ — * dS 



(85) 



(86) 



(87) 



4.3 Stability Criteria 

Tayler (1973) derives stability conditions for toroidal fields in a normal star in cylindrical coordinates using the energy principle 
given by equation (|78p . The equivalent conditions in spherical coordinates are given by Goossens & Veugelen (1978). We will 
now proceed to derive stability criteria for toroidal fields in a type II superconducting star, along the same lines. We will take 
the magnetic field to be given as a function of density and magnetic induction, H — H{p, B). This will allow us to consider 
both the strongly type II superconducting case and the normal case simultaneously. We will closely follow the notation of 
Goossens & Veugelen (1978) in order to facilitate comparisons. 

It is clearly sufficient for stability to show that the integrand of the energy of the perturbations is positive throughout 
the region of integration, 

SW=^ J £dV > if £>0 . (89) 

Even if £ becomes negative in a small region the system is unstable. Define £p and fmag as the integrands of SWp and (5Wmag, 
i.e. £ = £p + £mag. As in previous works (Bernstein et al. 1958; Tayler 1973; Goossens & Veugelen 1978; and Roberts 1981) 
we will drop the gravitational potential perturbation term in £p. The hydrostatic and magnetic parts of the energy are then 
given through equations (I74|) and (|83p . respectively, 



= 7P(V ■ if + a ■ Vp)(V ■ {) - (^ ■ V0)(V ■ Pi) 

1 r 1 (90) 

f mag ^ ^\SH -SB -SB -ix {V X H)^+ SpA^Pu 

The azimuthal angle ip does not explicitly appear in any of the coefficients in these equations, so that we can expand the 
components of the Lagrangian displacement as, 

= i?(r,e)e*'"^ , ^e = S(r,0)e™^ and = iT{r,6)e'"'^ . (91) 

Here m is an integer. Since only the real parts are significant, the scalar multiplications and vector dot products are to be 
treated as Z ■ Z* where Z* stands for complex conjugate. It will be of great notational convenience to define an operator A 
of a scalar argument u = u{r,8), 

A{u) = Rdru + . (92) 

r 

This is simply the directional derivative along the Lagrangian displacement, ^ • Vu = A(u)e""'^. We will find it convenient to 
redefine the tp component of the Lagrangian displacement as. 



Toroidal magnetic fields in type II superconducting neutron stars 17 



Also define, 



D = 



drjr'R) dgjSsme) 



+ 



r sin 9 



- T = Do - T 



(94) 



which is simply the divergence of the Lagrangian displacement, V ■ ^ = De"""^. Note that Do is independent of T. Using these 
definitions, we can express the hydrostatic part given by equation H90|) as. 



= 7pD2 + [a(p) - pA(<^)] D - A(p)A(0) . 
The equations of equilibrium for the unperturbed background state are given by equation (|8}, 

drP + pdr<f> = ~^dr {^^) " P^r 



dep + pde 



'^^"e sine — 
sm V vB 



pde 



dp 



(95) 



(96) 



Note the notational convention for partial derivatives that we will employ for the remainder of this section: derivatives with 
respect to coordinates x will be shortened as dx, while derivatives of the magnetic free energy F with respect to p and B 
will be explicitly written. Using these equations we can eliminate the pressure gradient in £p and rewrite it in terms of the 
gravitational and magnetic forces. Using the definition of the operator A from equation (I92p . we have, 



A(p) = -pA{(P) - pA 



OF 



BA 



dF 



B- 



dF f R + Scotl 



dp J \dB J dB \ r 

Next, consider the magnetic part of the integrand given by equation (|90|) . Using equation (|86p for 5H ■ SB, we have 



(97) 



47r 



SHSB+^ (\SBf - {5Bf 



SB ■ ^ X {V X H) 



- SpAipu 



(98) 



At/jii and SH are given through equations (|84p and (|87p . respectively. We can also express the magnetic field in terms of the 
free energy through equation H = AndF/dB. The various terms in fmag can be evaluated using the relations given in 
equation (|85l) . In particular, 



\SB\^ - {SBf _ m^{R' + S^ 



and 



SB ■ i X {V X H) 



B2 r^sin^e HB 

where we define the following auxiliary quantities. 



XY + TY , 



A{B) R + ScoH 



and 



B r 
The magnetic part can then be written as. 



Y = 



A{H) R + S cote 
H ^ r 



dF 

i^mag = B — 



m'{R^ + S'' 



XY -TY 



d'^F , d'^F,^^,2 d'^F ^ a^F, ^ 

^SpSB + —iSBf + ^SpAB + —SpAp 



(99) 



(100) 



(101) 



where 
SB 

"b 

AD 
D 

S_p_ 
P 

Ap 

P 



Do- 



R + Scote 



D 



Mp) 



(102) 



-De ^ 

We will next consider the m = and m 7^ cases separately. 



4.3.1 Them = Case 

In this case T = from equation (|93p and the total energy can be written as, using equations (|95p and poip for £p and fmag, 
respectively. 



(103) 
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where Do is defined in equation (|94|) . We have, in terms of the operator A defined by equation (|92|) . 



ICo = 'yp + B 
JCi = -2pA{(j)) - 2 
tC2 = -A(p)A(0) - 



dB 



dpdB 



R + S cote 



r 

^2 ! 



(104) 



A(5)|^ + 5A(B)_+BA(p)^ 



i? + S cot ( 



+ 



B 



dF 
"dB 



B 



2 d^F 
dB^ 



R + S cote 



All derivatives of R and S are included in Do- By completing the square we get, 



f = /Co Do + 



2K.0 



+ IC2 



4/Co 



(105) 



The first term is non-negative and the remaining terms form a quadratic in R and 5*, which is also the minimum value of £ 
with respect to Do, 



IC2 



4k;o 



= aoR^ + hoRS + coS^ 



(106) 



The subscripts in the coefficients stand for m = 0. Define the following auxiliary quantities. 



Ui 
U2 
Us 



-{B 

r 

1 



'dB 



28'^ F , „ d^F \ 



dpdBj 



dF d'^F d'^F \ 



dF 

4r ( deB— + BdeB 



dB 



dB- 



+ 



d^F \ 

dpdB ) 



(107) 



1 



dF 



B—+B 



,d^F 



dB dB^ 

We then find that the coefficients in the quadratic are given by. 



ao = 



bo 



Co = 



-drpdrCj) -U1+U3- -ji- I pdr(j) + Ul) 
/Co 



drpdg(f) dgpdrCj) 



Ui cot e-U2 + 2U3 cot e 



2 

/Co 



pdr(t) + Uo 



pdecp 



+ Uo cot t 



(108) 



depdefj) 



U2 cot e + Us cot^ G - — 



/Co V 



^P^'t + uocote'"' 



A sufficient condition for stability is that the quadratic form be always positive throughout the integration region. This 
corresponds to the following conditions, which are not all independent, 



a > , c> 



and 



b < 4ac 



(109) 



When these conditions are satisfied the star is stable, therefore these are sufficient conditions for stability. If we can show 
that the star is unstable as soon as one of these conditions is violated, then we will have shown that the conditions are also 
necessary for stability. For the m — case it can be shown that the interchange instability sets in when these conditions fail, 
as will be proven in a later section. Therefore, these conditions are necessary and sufficient conditions for the m = case. 
However, the same will not be true in general for the m 7^ will be discussed later. 

One way of deriving these conditions is to consider the minimum value of the quadratic form Q — aR^ + bRS + cS'^ with 
respect to S (or equivalently, R). For a minimum we need dQ/dS = and (fiQ/dS^ > 0. Substituting the value of S that 
minimizes Q and requiring that Q > we get the condition b^ < 4ac, while the second requirement gives c > 0. These two 
conditions then imply the third, a > 0. 

We can now consider special cases. In the strongly type II superconducting case the magnetic field is a function of density, 
H = H{p) and the magnetic free energy is given by equation ^ as F — HB/in. In particular, consider a power law of the 
form H (X p"^ . From equations (I104|l and (|107|) . we have. 



(UO) 
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so that the coefficients become, 

ao = -drp dr(t> h J - — pdr(l) H 

47rr 47rr^ K-o \ 47rr 

a.p^efli depdr(t> dr{HB) _ ^^(//B) , HB _ 
00= ^ cot — ^ h - — 7 cot 61 

CO = -^-^cot.+ i^4cot^.-4.fpa.^+i^±ll^cot.)^ 
47rr^ 47rr^ r^K-o \ 47r y 

On the other hand, in the normal conducting case the magnetic field and induction are equal H = B, and the free energy 
is F = B^ /Stv, so that from equations p04|) and (|107p . we have, 

^o = 7P+^, f/o = — , (71 = ^, C/. = ^ and C/s = ^ , (112) 
and the coefficients are given by, 

. . , BdrB B^ ^ f f,,.^ V 

drpde<t> dgpdr(i> BdrB BdgB B^ 

oo = cot — — H cot 

r r Znr 2nr^ nr^ 



2 / o . ^2 W „ , B^ 



(113) 



CO = 5 — 5- cot 6'+- — -cot ^71^ \ pOe(p + — cot t 



These are the same as the results given by Goossens & Veugelen (1978)0 



4.3.2 Them^O Case 



When m 7^ 0, the hydrostatic and magnetic parts of the energy are given by equations (|95|l and (|10ip . respectively. In this 
case, the integrand £ = £p + fmag is quadratic in the rescaled ip component of the Lagrangian displacement T, defined by 
equation (|93[) . and does not contain any derivatives of it. Therefore, we can write the energy as. 



-9 - OF 
£ ^ £0 + aT^ + f3T + B^ 
oB 



where £0 is the energy for the m = case, given by equation (|103p . and we define, 



(114) 



" = 7P + P 



/3 = -2 7p + 



„ d^F ^d^F 



Do + 2pB 



d'^F fR + S cote 



dpdB 



(115) 



2pA{^) 



£0 is independent of T. We therefore have S£ jdT'^ = 2a. The 7p term in a will be the dominant term for the cases of interest 
to us, so that d?£/dT^ > 0, and consequently £ can be minimized with respect to T. Setting d£/dT = we get the value 
that minimizes the energy, T = -~(3/2a. Substituting this back into the energy we find the minimum as. 



£^£0 



^+B^ 
4a dB 



(116) 



As was done in equation (I103|l for the m = case, we can once again group together terms of different order in Do, defined 
by equation lf94)) . 



£ = CoDo' + £iDo + £2 



(117) 



^ Note that there is a typo in equation (13) of Goossens & Veugelen (1978). 
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For notational convenience, define a set of auxiliary quantities, 
Vo ^ a ' 7p + P-B „ a D + P 



pae<^ pB cote d'^F \ 



dpdB J 



and. 



W2 = - Uo cot e» + v;)V2 

r 

where a is defined in equation (|115|l . and Uo is defined in equation (|107|l . Also invoking the definitions of Jd from equation 
(|104p . we have, 



/:o = /c„-Vo^ = B^^--^pB^^^^^ 

£1 = /Ci + 21/o(Vi7? + V2S) = 2(Wii? + W'zS') (120) 

C,=JC2-{V^R + V2Sf + B^ + 

oB 

Rearranging the terms we get. 

Note that £0 is not necessarily positive, so unlike in the m = case, it is not obvious that the first term is positive definite. 
In fact, for the strongly type II case where the free energy is of the form F — H(p)B /in, we have jCq < 0. On the other 
hand, for the normal case F — B^/Svr, so that Co > 0. For negative Co the system is unstable since we can find displacement 
fields with sufficiently large derivatives which will make the Do term dominant in the energy. Therefore, for stability we must 
require £0 > 0, or using the definitions of equation (|120[) . 



d^F 



■V 



dB-^>y''d^) /[^^^pw' 



This is a necessary but not sufficient condition for stability. This is related to what we will refer to as the Muzikar-Pethick- 
Roberts (MPR) instability (Muzikar & Pethick 1981; Roberts 1981), which we will discuss in more detail in a later section. 

Another way of looking at equation (I121|l is to say that when Co > 0, the energy can be minimized with respect to Do- 
The minimum is a quadratic in R and S, just like equation (|106p for the m = case. 

The coefficients are given as, using the definitions of Ui, Vi and Wi made in equations (|lU7p . (|118|l and (11191) 



2 

= a^ie + br^RS + c^S' . (123) 



flm = —drpdr4> — Ul + U'i — Vl 



2 rn^B dF 
^ r2 sin^ e dB ~ £^ 



drpd0(t> dgpdr(p , „ rr , r,rr .n OT/ 1/ '2WlW2 

(7i cot 6 —U2 + 2(73 cot 9 — 2V1V2 (124) 

r r Co 

depdocp 2 2 , rn^B OF Wi 
^ U2 cot 9 + U3 cot 6 — V2 



r2 sin^ e dB Co 

This quadratic is positive if the coefficients satisfy the conditions listed in equation pU9p . However, the system will be definitely 
stable only when Co > 0. On the other hand, if these conditions are violated, i.e. if the quadratic is negative, then the system 
is unstable regardless of the sign of Co- Also note that, clearly, the |m| = 1 case is the worst instability, as noted previously for 
the normal case by Tayler (1973) and Goossens & Veugelen (1978). On the other hand, when £0 < the energy is maximized 
with respect to Do, and it is always possible to find a Lagrangian displacement field with sufficiently large derivatives that 
will make the system unstable. 

The coefficients for the strongly type II case can be obtained by setting F = HB/An. On the other hand, for the normal 
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case we have F — B^/Svr, and the coefficients reduce to, 

am = -drpOr(p " T + 



7P 27rr 27rr^ 47rr^ 

drpde(l> depdr(l> 2p^dr4>de(t) BdrB BdeB B 



cote-— ^ J cote (125) 



depde<t> {pdB<i>? BdeB _ B^ , m^B^ 

■ cot B — - — - cot + 



J.2 jpr'^ 2iTr^ 27vr^ inr^ sin 6 

These are the same as the resuhs given by Goossens & Veugelen (1978). 

In the next two sections we wiU consider the special cases of the completely normal conducting star and the strongly type 
II superconducting star with H oc p. The coefRcients a, b and c (given by equations 11081 and 1 124[) have hydrostatic terms that 
are of the form drP and dep, and magnetic terms of the order of the magnetic free energy F. The radial dependence of the 
background quantities arises from the much stronger hydrostatic forces, while the dependence arises as a result of magnetic 
forces. Therefore, dep ^ F <^ drP- We will calculate the coefRcients to first order in the magnetic energy, which is much 
smaller than the hydrostatic terms. We will assume that the perturbations and the background state have the same index, 
thus neglecting buoyancy effects. If we include buoyancy, then to leading order, the coefficient a will be a buoyant term, c 
will be a purely magnetic term, and b will be the product of a buoyant term and a magnetic term. Thus, b^ <^ Aac, and the 
stability conditions (given by equation I109|l will reduce to a > and c > 0. The first condition is necessary for stability to 
buoyancy, and the second is the same condition on the magnetic field as without buoyancy. We will consider the effects of 
multi-fluid composition in more detail in future work. 



4.4 Stability Criteria for a Normal Star 



We will now examine the stability of a particular magnetic field configuration in a normal star. The equilibrium equations in 
this case are, from equation (|96[l. 



drP + pdr4> ~ — 
dep + pde4> = - 



BdrjBr) 

4:Tvr 
B de{B sine) 

in sin 9 



(126) 



Let po, po and (po refer to the hydrostatic equilibrium in the absence of magnetic fields. This equilibrium is spherically 
symmetric and is simply given through. 



drPo + Podr<t>o — 



(127) 



The difference between po, po and (po and the corresponding quantities p, p and (j> in the presence of magnetic fields is of the 
order of the magnetic pressure ~ B^, which we assume to be small compared to the hydrostatic pressure. Therefore, using 
the equations of equilibrium we can expand equation (|113|l for m = to lowest order in B^, 



fflo 



bo 



Co 



" 3B^ 



BdeB 



din po 
dlnr 

cote - 



+ 
BdeB 



3-B-' 



BdrB 



din po 
dlnr 



cote 



B^ 



d In po BdrB 



■ cot 



din r 

2 o 



2Tvr 



cot ( 



BdrB B^ 
2nr 

BdeB 



2-Kr^ 



+ 



2-Kr'^ 
B' 



■ cot f 



(128) 



2-Kr'^ 2-Kr^ 
On the other hand, for m = 1, we have, from equation (|125p 



- 2 + 



BdeB 
' 2nr^ 



din Po 
dlnr 

din po 
dlnr 

cote - 



B^ BdrB 

+ —r— 



B^ 
B^ 



+ 



B' 



Atvt'^ sin e 



cote + 



27rr2 



■ cot 



BdeB 
47rr2 

B^ 



BdrB 
2Tvr 



cote 



2Tvr^ 



■ cote 



(129) 



471-f 2 



sm 



We will now consider a specific example. Let the equation of state be given by a 7 = 2 polytrope, where the background 
density profile is p = pc sin x/x, in terms of the dimensionless radial coordinate x = r/ro- Assume a magnetic field of the form 
given by equation (|4T|) . 



B{r,e)^Bo 



(n+2)/4 



i/2 



T./2 . 



= BoX^ 



-2)/4 . (n+2)/4 . n/2 , 

" snv " xs\n ' ( 



(130) 
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where n ^ 1. Then, for m = 0, the coefficients become, from equation (|128p , 
16nri 

bo « -^{2-n){l + xcotx)x^"''^^/''sm'-"+^^^'' xsin"-^ 9cose (131) 

-"o /o \ (n-6)/2 • (n+2)/2 • n — 2a 2 a 

Co ~ ;r(2--n)x^ " sm^ " a;sm &cos 

Note that 60^ = 4aoco, so that the quadratic forms a complete square, i.e. aoR^ + boRS + coS'^ — ao{R + boS/2ao)'^ ■ However, 
for n > 2, we have ao < and co < 0, and the conditions for stabihty (equation I109[) are violated. Thus, only fields with 
1 ^ n ^ 2 are marginally stable for m — 0. 

On the other hand, for m = 1, we have, from equation (|129|l . 

^ [4 - (n + 2) (1 + a; cot sin^ ^1 x^"-^)/^ sm^"+^^^^ x sin"-^ 6 

^" ' +2)(l+:Ecota;)a;'"-'')''"sin'"+"'/''a;sin"~'ecos6' (132) 



47rr^ 

Bo r / I ri\ 2 fli (ri — 6)/2 . (n + 2)/2 • n-2 , 

— J [1 — (n + 2) COS 0\ x^ " sm^ " xsm I 



Since a™ and Cm become negative in some regions, they violate the stability conditions given by equation (|109|l . Consequently, 
the normal magnetic field is unstable for m — 1. Thus, we might expect n = 1 models with both normal and superconducting 
regions to be unstable. Poloidal fields may stabilize the star, as in normal conductors (Tayler 1973; Wright 1973; Braithwaite 
& Nordlund 2006), and we consider adding them in a following section. 



4.5 Stability Criteria for a Superconducting Star with H oc p 



We will now consider the strongly type II superconducting case with H oc p (i.e. a = 1) in more detail. In this case F = HB/in, 
and the equations of equilibrium (equation [96]) explicitly give. 



drP + pdr(f> - 



Bdr{Hr) HdrB 



^ . , HB HdeB 



(133) 



Using these equations as well as the equation of equilibrium in the absence of magnetic fields (equation [127]), we can expand 
the coefficients for m = (equation lllip to lowest order in HB, 



ao 



bo 



Co 



HB ( dlnpo 
d\nr 



27rr2 
?,HB 



+ 



3HB dr{HB) 



cote - 



deiHB) 



4j^j.2 

deiHB) 



cotO + 



HB 



dlnpo dr(HB) 



dlnpo drjHB) HB 
dlnr 47rr Aixr"^ 

de{HB) HB 



dlnr 



cot f 



47rr2 



+ 



27rr-2 



cot ( 



(134) 



47rr^ 47rj-2 

For a 7 = 2 polytrope we have p — pcSinx/x. Consider a magnetic field of the form given by equation p5|l . for ct — 1, 



B{r,9) = Bj-^] — sin" 6» = Bo sin" a: sin" e 



where n J5 1. In particular, we get, from equation (|134|) . 
Co 



HcBo ,^ , _3 . n + 1 

-(1 — n)x sm - 



47rr-; 



(135) 



(136) 



For all n > 1 this is negative, thus immediately violating one of the conditions for stability (equation [T09|. For n = 1 all three 
coefficients vanish to lowest order in HB, implying that the magnetic field is marginally stable. In Appendix A, we show that 
this result is true for any //(p, B). 

For m 7^ 0, we have f equation 1 120|) . 

2 



do 



1 



(HB\ 
7P \ 47r ) 



< 



(137) 



which implies that even if the conditions given in equation p09p are met the system will still be unstable. This is the MPR 
instability and will be discussed in a following section in more detail. 
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4.6 Interchange Instability 



In this section we will show that the m = stability conditions correspond to the stability criteria for the interchange 
of two magnetic flux tubes, as demonstrated for the normal case by Tayler (1973). Consider two axisymmetric flux tubes 
located at coordinates r, 6 and at r + Sr, 9 + 59, and having volumes V and V + 5V and corresponding cross-sections A and 
A + &A, respectively. We will assume that the interchange is adiabatic so that the mass pV , magnetic flux BA and pV^ are 
all conserved. 

Let the pressure, density and magnetic induction of the two tubes initially be, 
at r, ^: p p B 

(138) 

a.t r + 6r, 9 + 59: p + Sp p + 5p B + SB 
After the interchange the corresponding quantities are, defining a cylindrical radius by w = r sin 9, 



at r, 9: 



at r + 5r, 9 + ( 



{p + Sp){V + 5V)'' {p + Sp)(y + 5V) {B + SB){V + 6V)uj 



pV 



V 
pV 



V{-co + 5vj) 
BV{-w + 5zo) 



(139) 



{y + 5V)i V + 5V {V + 6V)-cu 

The total energy is the sum of internal, magnetic and gravitational energies. Without loss of generality, we can take the 
zero of the gravitational potential to be at r, 9. Prior to the interchange, the energy is, 

pV (p + 5p)(V + 5V) 



+ F(p, B)V + F{p + Sp,B + 5B){V + SV) + {p + 5p){V + SV)5<j} 



7 — 1 7 — 1 

Here F is the magnetic free energy. After the interchange, we have 

pV 



{p + 5p)iV + 5V)-' , 

Ef = 7 -.xw. -, 1- 



- + F{pi,Bi)V + F(p2, B2){V + 5V) + pV5(l> 



(140) 



(141) 



(7-l)VT-i (7- 1)(V + 5V)^- 

Here p\ and B\ are the new density and induction at r, 9, and p2 and B2 are the corresponding quantities at r -\- 5r, 9 + 59 
(equation [139]). We need to calculate the energy difference resulting from the interchange to second order, 

AE = Ef - -E, = AEp + AEm ■ (142) 

Here for notational convenience we denote by AEp the change in the internal and gravitational energies, and AEm is the 
change in the magnetic energy. To second order we have, 



^Ep « JP^^^y- + {Sp - p5(f))5V - V5p5(f, . 

Using the equations of equilibrium (equation [96]) we have. 



5p - 



dF 5vj 



d^F 



d^F 



d^F. 



The magnetic term in the energy change is lengthy. First, note that, 

d^F 



(143) 



(144) 



(145) 



Fip + 5p,B + 5B) . Fip, B) + -5p + -5B + - — (5pf . ^ 

We can write the magnetic terms in Ef (equation 11411) as F{pi,Bi) = F{p + 5pi,B + 5Bi), so they can be expanded in a 
similar fashion. Here we have, to second order. 



{p + 5p){V + 5V) 
Pi- T7 - P = P 



5pi 

5p2 — p2 — p — 
5Bi = Bi-B = 



V 



pV 



V + 5V 



5v; 

' V 



5p 5\^ 5p5V_ 



5V 



V 



{B + 5B)(y + 5V)txj 
V{vj + 5vj) 



B 
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Using these and equations (|143p and (|144p we can write the energy change as. 
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where, 

Since A^o > for the cases of interest, the change in energy can be minimized with respect to 5V/V . The minimum of the 
energy becomes, 

^^M.-^. (149) 

The small quantities need to be expanded only to first order, 
Szu — S{r sin 6) = Sr sin 6 + rS6 cos 9 

(150) 

Sp = Srdrp + Sddgp 
and similarly for B and (p. The energy can then be written as, 
A F 

— ^ ao{Srf + borSrSe + cor^SOf . (151) 



ao, bo and cq are the same as in equation (|108|) and the conditions for stability are the same as in equation (|109|l . In fact, 
the same conclusion could have been drawn by comparing equation (|148p to (|104p . Thus, we have shown that the m = 
stability conditions are the same as the conditions for stability under the interchange of magnetic flux tubes. In other words, 
the interchange is the worst instability for m = 0. 

4.7 The Muzikar-Pethick-Roberts (MPR) Instability 

In this section we will derive the criteria for the instability discussed by Muzikar & Pethick (1981) and Roberts (1981). Using 
equation ((2)1 we can write the magnetic stress tensor as (equation [l|, 

= (F - pF,p - BF,b) + BF,bB,B, . (152) 

Consider perturbations around a state of uniform density p and uniform magnetic field B — Bz. The Lagrangian displacement 
associated with these perturbations is, 

4(r, t) = ^ exp(ife • r — iut) . (153) 
In this case, we have, 

= V X (I X B) = iB{k^i -k-iz) 
SB = B ■ SB = -iBk± ■ 

(154) 

SB = B'\SB - SBB) = ik,i^ 

Sp = -V- (pi) = -ipk ■ $ 

Here _L means perpendicular to z. 

The magnetic force density is, from equation p52|) . 



= V • o- = - ipF,pp + BF,pb) Vp - {pF,pB + BF,bb)VB + B ■ V{F,bB) . (155) 
Since the background quantities are constant the perturbation in the magnetic force becomes, 
Sf^ = -(pV,pp + pBF,ps)fc(fe-|)- (pBFpB + BV,BB)fe(fcx-$x) 

(156) 

+zk, [pBF,pB{k ■ + B^F,BB{ki_ ■ |^)] - BF,Bkk^ 



Since the background state is symmetric with respect to z we can choose k — zkz + xk^. With this choice equation (I154|l 
becomes, 

SB = iB{kz^xX + kz^yij - kx^xz) 

SB = —iBkxix 

(157) 

SB = ikz{^xX + ^y-y) 
Sp = —ip{kx^x + kz^z) 
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The components of the magnetic force become, 
((5/,n)^ = [kl {p^'Fpp + 2pBFpB + B^Fbb) + klBFe] - ^.k,k, {p^F,pp + pEF^ps) 

{SU)y = -iykiBF,B (158) 
iSfm)z = —^xkxkz {p^F,pp + pBF,pb) — izklp^F^pp 
In addition, there is a pressure restoring force, Sfp — — V5p = —^pk{k ■ ^), or in components, 

{Sfp)x = -'yp{kUx + kxk^z) 

[Sfp)y = (159) 

{5fp)z = -')p{kxkz^x + kl^z) 
We will neglect gravitational forces, so that the equations for the perturbations become, 

- puj^i = Sfp + Sf^ . (160) 
From equations (|158[) and p59[l it follows that the equation for completely decouples from the equations for ^j, and ^z, 

Pi^^^y = klBF,B^y . (161) 

This implies that one pair of modes has S,x ~ ~ and 7^ with lo'^ = kj.BF,B / p- These modes are the generalization of 
the Alfven modes. The remaining modes are given through, 

puj^ix = ix \kl hp + p'F,pp + 2pBF,pB + B^F,bb) + k^BF.B] + izkxkz hp + p^F,pp + pBF,pB) 

2 22 (^^2) 

pijj'^£,z = ixkxkz (7P + p F,pp + pBF,pb) + £,zk^ [-yp + p F,pp) 

From these two equations we get the characteristic equation for the modes, after some rearrangement, 

p^oj^ — pLo^Eo + f 1 = , (163) 
where, defining k^ — k^ +kl, 

jBF.nR + B^Fhh) + kl (BFr + p'^F.nn) 

(164) 



£0 = k^JP + kl (p^F,pp + 2pBF,pB + B^F,bb) + kl (BF,b + p^F,, 



£1 = klkl {jpB'^F^BB + p'B''F,ppF,BB - p^B'^fIb) + ktBF,B (tp + P^F.pp) 
In the absence of magnetic fields, we have, defining jp — pel, 

p^Lu"^ — p'^uP'k^cl — , (165) 

which has two roots: a;^ = and ui^ — k^cl. The latter corresponds to sound waves. In the cases of interest, the magnetic 
terms will be much smaller in comparison to the pressure terms, so that one of the roots will have cj^ ~ k^cl and therefore 
will be definitely positive. Since £1 is the product of the two roots, the condition for stability is £1 > 0, which for kz ^ 
becomes, 

kl {-ypB^F^BB + p'B'^F.ppF.BB - P^B^F^b) + klBF.B (7P + P^-F.pp) > . (166) 
For sufficiently large kx, or more precisely when klBF^BB S> klF^B, this reduces to, 

^,BB > ; — 5-= — ~ . (167) 

IP + P Fpp IP 



This is exactly the same condition for stability as in equation (I122p . When pressure dominates, it is also of the same form as 
the condition given by Roberts (1981). From equation (|163|) it follows that the potentially unstable modes are given through, 

2^ £1 klkl f . P^B^FIb \ , fcto^ 

~ T2 — ~ ~TT' ^ ^'BB — + T^BF^B ■ 168) 

/e^7P k'' \ IP J k-^ 

The magnetic free energy in the strongly type II case {H ^ B) can be written as (Tinkham 1975; Muzikar & Pethick 
1981), 

^=^^ + \/ = tI T exp . (169) 



47r V 327r3 \a J A 

Here <I>o = hc/2e is the fiux quantum (n$ = B/^o is the fiux line density per unit area), A = {nipC'^ /A-KUpe'^y^^ is the London 
penetration depth, Up is the number density of protons, and a is the distance between flux lines in a triangular lattice. 



"= 37 \B 



(170) 
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The magnetic field strength in this case is (Tinkham 1975; Easson & Pethick 1977), 

where ^ = h^kp/i^rnpA is the coherence length, ^ <^ X; kp = {Sn^np)^^^ is the Fermi wave number of protons, and A is 
the superconducting energy gap. The first term in equation (|169[l is the energy of an isolated flux line, and the second term 
arises due to the interaction between flux lines. Note that only a depends on B and only the interaction term contributes to 
F^BB- Also note that cx 1/p when the proton number density is proportional to the baryon number density, as suggested 
by Baym et al. (1971). Defining a new variable by u = a/ A we have (equation [169}, 



4, +E{p)u-^^'e- where iJ(p) = ^^^. (172) 
Then, introducing an auxiliary function f{u), 



B^Fbb = ^(u-'^^+2u''/^ + 



^u-'/'^ e- = E(p)i{u) . (173) 

Only the first term needs to be retained when u ^ 1, i.e. when the spacing between fiux lines is large compared to the 
penetration depth. In the same limit, we can also approximate, 

pBdH gHB , d\nH 

pBF^pB ^ — = — where a = — . 174 

47r dp An d\np 

Using these equations, we can write the condition for instability as, from equation ([TSTJ, 

u'f{u)<J-^^. (175) 

Note than when a — 0, i.e. when H is independent of p, there is no instability. Thus, it does not arise in a normal medium. 
Moreover, cr > is not required in order to have an instability, contrary to the conclusions of Muzikar & Pethick (1981). 

We take the magnetic field strength to be // ~ lO'^^ G, and the typical density in the superconducting region to be 
p ~ 3 X 10^* g/cm^, which corresponds to a pressure p ~ 4 x 10"^"^ erg/cm^, for a 7 = 2 polytrope and a radius Ri, ~ 10 km. 
We also take a = 1. From equation (|175|) it follows that instabilities arise for u > Uo where Uo — 20. Using equation (|170|l 
and the definitions of A and $0, we can find the largest magnetic induction which is unstable, 

B < = 1.15 X 10^3 f^^V-Vc, (176) 

The proton number density Up is a function of the baryon number density, and for nj, ~ 0.2 fm~"^, we have Up ~ 0.01 fm~^ 
(Elgar0y et al. 1996; Zuo et al. 2004). 

For toroidal fields z is along the </> direction, so that for modes we have exp{ikzz) — exp(im(;/>). We can take kz ~ m/Ri, 
for a star of radius R^ . The condition given in equation p67p can lead to instabilities when the perpendicular wave vector kx 
is sufficiently larger than kz- Using equations (|172|l and (|173[) . we get, 

kl BFbb A-KE{p)f{u) _ 3V2^u\f{u) ^ a^HB 

kl F,b ^ H{p)B ln(A/e) ~ 47r7p ' ^ ' 

where the last inequality follows from the condition for instability f equation 1 175|) . The length scale of the instabilities is small 
compared to the size of the star; for a 7 = 2 polytrope. 



= k^ <^ — t — « 3.1 X 10 cm . (178) 

m y 47r7p mpn 

Here H15 = H/IO^^G, B12 = B/IO^^G, and pi4 = p/10^*g/cm^. From equation (I168|l we can estimate the growth rate of 
the instability, using yp — pc^, 



kzBF,pB 



m\a\HB 



AnpCsR* 

The corresponding growth timescale is. 



(179) 



1^3.7x10^ s. (180) 

ui m\a\Hi5Bi2 

Here Ra = Ri,/10^' cm. Note that m = is stable. The unstable modes will be dissipated if the kinematic viscosity of the fluid 
is, 

^>|-26i^cmVs. (181) 
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This value is well below the estimated values of the viscosity in a neutron star, which are typically in the range rj ~ 10*"^ cm^/s 
(for a review see Andersson, Comer & Glampedakis 2005). 

Note that similar results will hold for poloidal fields, except that in this case fcz > 1/_R^ will depend on both the number 
of radial nodes and the angular momentum quantum number of the mode. Simple linear analysis along the lines outlined by 
Hide (1971) reveals that the growth rate of the MPR mode will not be strongly affected by buoyancy, but the condition for 
stability will be modified. Moreover, due to the local nature of the mode, it is likely to be unaffected by rotation. 



5 NEARLY TOROIDAL FIELDS 

In normal conducting stars, the presence of poloidal components in addition to toroidal components may help stabilize the 
magnetic fields (Tayler 1973; Wright 1973), which has also been confirmed by recent numerical simulations (Braithwaite & 
Nordlund 2006). Moreover, pulsar observations reveal the presence of a dipole-like field in the neutron star magnetosphere, 
implying that a poloidal component of the magnetic field must exist. The treatment of fully poloidal fields is considerably 
more complicated and will be discussed in a subsequent paper. The complication arises as a result of the fact that in the 
poloidal case the direction of the magnetic field is not known, and must be computed numerically (Roberts 1981). 

In this section, we will consider the case when there is a small poloidal component in addition to the much larger toroidal 
field. We will evaluate the constraints on the shape of the poloidal component that result from the restrictions that the 
magnetic force per unit mass be expressible as a gradient of a potential and that V ■ i? = 0. We will then consider the 
boundary conditions that must also be satisfied. Let us assume that the direction of the field is given by, 

n^4> + e, (182) 

where e is a poloidal vector and |e| ^ 1. In what follows, we will retain only the first order terms in |e|. 

The form of the magnetic field inside the superconductor is H — H{r,6){(j) + e) and the current density can be written 
as the sum of toroidal and poloidal components, so instead of equation (|13p . we now have, 

J — J tor J pol 

47rJtor „ „i V(i/r-sin6l) X 

= V X H4> = — ^ (183) 

An J pol _ rr 

— = V X He 

c 

Note that Jtor (due to the toroidal magnetic field) is a poloidal field and Jpoi (due to the poloidal magnetic field) is a toroidal 
field, i.e. Jtor -L 4> and Jpoi || 4>. Taking the induction to be i? = B{r,8){(j) + e), the first term in the force density, given by 
equation (0, becomes, 

JxB Jtoi X B4> JpoixB4> Jtor X Be 

= \ 1 . (,1841 

c c c c 

The first term is due to the toroidal field, and the second and third term are due to the presence of the small poloidal 
component. Since Jpd is a toroidal field the second term vanishes. On the other hand, the third term is a cross-product of two 
poloidal fields, and therefore is a toroidal field. However, we require the toroidal force density to be zero, so it must vanish. 
This means that e || Jtor, or equivalently, in terms of an arbitrary function A, 

B{r,e)e = \{r,e)Jtor ■ (185) 

Thus, the force is of the same form as in the purely toroidal case, and in order for it to be a gradient, the induction B must 
still be of the form given by equation psp . We get a condition on the unknown function A from V ■ B = Jtor • VA = 0, 

47rJtor-VA ^-VXxViHrsme) ^ 

c r sm 6' 

This equation is satisfied by functions of the form, 

X{r,e) = A(_ffr-sine») . (187) 

Thus, the poloidal vector e is given by equation (|185p . using equation (I15p for B and equation (|183p for Jtor, 

A (ffr sin 61) Jtor H^-,-.^^ . ■> 

^ = 1 . auu ■ Q\ = VA Hrsm^ x . 188 

4:Tv pr sm J [HrsmOj p 

In a normal conductor, we have, setting H = B and using equation (|18|) for B, 

fl{Br sin 6) J tor r-,~/r, ■ n\ 1 

e = — = V^{Brsm6) x </> . (189) 

B 

Here fj, and jl are arbitrary functions. 
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5.1 Boundary Conditions 

Neglecting second order terms in the small quantity |ej in the magnetic stress tensors for the normal and superconducting 
regions (equations |3] and |4|, the boundary conditions for the continuity of stress (equation [25]) become, 

— Sps + CTrr.s = —Spn + 0"rr,n and (Jrifi^s ~ (^r<li,n ■ (190) 

The rr components of the magnetic stress tensors are the same as in the purely toroidal case fequation l29|l . so the first equation 
is the same as before (equation [27} . However, we now have the second equation, which explicitly gives, using equations ^ 
and ^ for the stress tensors, 

(4>-H){f-Bs)^{4>-B^){f-B„). (191) 

We also have the additional boundary condition on the continuity of the normal component of the poloidal magnetic induction, 
which follows from Maxwell's equations, 

f Bs=f B„ . (192) 
The last two equations imply that we must have, 

$ ■ H = 4> - Bn i.e. H = Bn . (193) 

This is equivalent to the requirement for the continuity of the </> component of the magnetic field in the absence of surface 
currents (equation I20p . However, as was previously discussed, this is inconsistent with our assumption that _ff is a function 
of radius up to the boundaries of the superconductor. This assumption now requires the presence of a discontinuity in the 4> 
component of the magnetic force, although the forces within the superconducting and normal regions have no such components. 
This is an artifact of the incomplete description of the transition boundary, which we have treated as discontinuous. A more 
realistic treatment should impose zero toroidal force everywhere. 

Incidentally, note that we cannot simply assume that the radial components of the poloidal vectors vanish at the boundary, 
which would also satisfy the above equations (equations ll91l and ll92p . This would imply that the functions A and p. in equations 
(|188p and (|189|) are constants, which in turn would cause the poloidal vectors to vanish everywhere within the normal and 
superconducting regions. 



6 CONCLUSION 

Our main goal in this paper has been to compute the distortion of a neutron star due to a toroidal magnetic field in its interior, 
assuming that the star is either partly or entirely a type II superconductor. Previous authors have estimated the order of 
magnitude of this distortion (Jones 1975; Easson & Pethick 1977; Cutler 2002), finding that it is enhanced by a factor H/B 
for given magnetic induction B and magnetic field H compared with the normal case (where H — B). In the strongly type 
II regime, H ~ 10^^ G, so that H/B ~ 10^/-Bi2 (Jones 1975; Easson & Pethick 1977). Such large enhancements could result 
in magnetic distortions e ~ 10"'' — 10~*, which are large enough to be important for neutron star precession (Wasserman 
2003) and possibly for gravitational radiation emission (Cutler 2002). These earlier works did not compute the structure of 
the magnetic field in detail. 

Here, we have paid closer attention to the requirements of hydrostatic balance and stability. The assumption of a barotropic 
equation of state, p = p{p), which ought to apply to a cold neutron star, severely constrains the structure of the toroidal field. 
Similar restrictions have been known for a long time for normal conductors (e.g. Prendergast 1956; Monaghan 1965). The 
restrictions arise because the magnetic acceleration must be a total gradient in hydrostatic balance. Under the assumption 
that the magnetic free energy F is a function of (matter or baryon) density p and magnetic induction B, we find that, for 
toroidal fields, we must require (equation I15|) . 

B{r, 6) oc pr sin ef{Hr sin 6) , (194) 

where / is an arbitrary function. Given this function, and F{p,B), we can compute H{p,B) — AndF/dB (equation [2]). 
Equation (|194p is then an implicit equation that can be used to find B(r,6) (assuming axisymmetry) . Similar constraints can 
be derived for poloidal magnetic fields, but are more complicated since the field direction must be solved for (e.g. Roberts 
1981 for superconducting, uniform density stars; we will consider superconducting, barotropic stars in a future paper). 

Our calculations have concentrated on neutron stars with a strongly type H regime where H is independent of B; our 
models allow for as many as two normal regimes interior or exterior to the superconductor. The main result of these calculations 
is equation (|7ip for the magnetic distortion, 

with (j>2{R*) /"^o ~ —2 in all cases, as is summarized in table[T] These results were computed for an equation of state p — np^ 
and H cx p (e.g. Easson & Pethick 1977). Calculations can be done in a similar way for other p{p) and H{p,B). 

Although we have separated the star into strongly type II and normal sectors for computing the deformations due to 
a toroidal field, we have noted that this assumption, while mathematically well defined, leads to sudden jumps in density 
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and magnetic induction at the boundaries of the superconductor. In effect, we have assumed that the magnetic free energy 
changes discontinuously from F — H{p)B /An in the type II superconductor to F = /Stt in the normal conductor. However, 
our formalism can be applied more generally to F{p, B) that varies smoothly from type II to normal, probably with inter- 
mediate domains of type I superconductivity. Such models ought to be free of discontinuities in p and B, but will still have 
rapid vaxiations in radially thin domains. In particular, we expect magnetic stresses to be approximately continuous across 
boundaries, so the magnetic induction _B„ in the normal regions will be larger than the induction Bs in the superconductor, 
Bn ~ [HBsY^"^ > Bs. Strong toroidal fields B„ ~ 10^^ '^ G (corresponding to // ~ lO^'^ G and Bs ~ 10^^ G) are needed for 
large distortions; toroidal fields _B„ ~ lO^'^ G imply Bs ~ 10^ G and therefore will lead to e ~ 10~^^. We have postponed 
considering models with realistic F{p,B), which would be more intricate mathematically, to later work. 

A toroidal field can be produced as a result of the winding up of the magnetic field early in the history of a neutron star 
(Thompson & Duncan 2001). The resulting field could be stronger than lO^'^ G. When the star has cooled down sufficiently, 
the superconducting shell forms. This would produce a large stress within the superconductor and the star would become 
dynamically unstable. This, in turn, would lead to a lowering of the induction inside the superconductor until stability can be 
restored. In equilibrium, the stresses within the superconductor and the normal regions will be comparable. In other words, 
the amplitude of the magnetic stress may be fixed by the original amplification of the toroidal field. The superconductor 
adjusts to the requirement of approximately continuous stress by lowering Bs- In this sense, the superconductor doesn't really 
amplify the stress. 

Magnetic fields not only need to be in magnetohydrostatic equilibrium, but they must also be stable with respect to 
perturbations. We have derived stability criteria from an energy principle for generic F{p,B). This is more general than 
the treatment of Roberts (1981), who assumed H oa p, and it also includes the normal case treated previously by Tayler 
(1973) as the special case H = B. In a completely type II superconducting star with H (x p and B oc sin" 9 (equation 
I35p . we find that only n = 1 is stable to m = (axisymmetric) perturbations. In fact, as we show in Appendix A, this 
is true for any magnetic field of the form H{p,B). For m 7^ all field configurations in a type II star are prone to the 
Muzikar-Pethick-Roberts (MPR) instability, found by Muzikar & Pethick (1981) and Roberts (1981), when B < 10^^ G. 
There is also a minimum wave number for instability, and it is very large: the MPR instability is a small scale instability. 
From a linear perturbation analysis around a uniform background, we find that the instability has a length scale ~ 10"'' i?*, 
where Ri, is the stellar radius, and a timescale ~ 10^ s. This timescale is relatively long compared to an Alfven crossing time 
tA = Ri,{'^-Kp/ H BY^^ ~ 3.5 Reipis / Hi5Bi2y^^ s, but short compared to a typical precession period of the order of a year. We 
have also argued that the MPR instability cannot occur for m = in toroidal fields: our linear analysis implies zero growth 
rate for modes with wave vectors entirely orthogonal to the unperturbed magnetic field. Because of the large wave numbers 
required for the instability, viscous effects, which cannot be studied via stability analyses from energy principles, could prevent 
it from occurring altogether. Our estimate is that a kinematic viscosity of ~ 10 — lOOcm^/s would be enough to shut off the 
instability; this value is smaller than most estimates of the kinematic viscosity in neutron star matter (Andersson et al. 2005). 

We find that normal toroidal fields with B oc sin"'^^ 9 (equation 141 1) are unstable for m = 1. Therefore, toroidal fields in a 
star with normal and superconducting regions will be unstable. Poloidal fields may help stabilize the stellar magnetic field, as 
has been found for normal conductors (e.g. Tayler 1973; Wright 1973; Braithwaite & Nordlund 2006). Moreover, the emission 
from radio pulsars additionally requires exterior, poloidal fields. Consequently, we have also considered nearly toroidal fields 
in which the field direction is t^ + e, where e _L <^ and jej ^ 1. Here, too, the form of e is not completely arbitrary: to maintain 
hydrostatic balance and eliminate toroidal forces, we find the requirement (equation [188]), 

tt2 

e = VX(Hr sin 9) x 4> , (196) 

P 

where A is an arbitrary function. We derived equation (|196p for type II regimes, but it holds elsewhere (in particular, in normal 
regions). We have seen, though, that when we assume discontinuous transitions in the magnetic free energy between type II 
and normal regions, there are discontinuities in the r(j> component of the magnetic stress tensor, implying a surface toroidal 
force. A more complete treatment with continuously varying F{p, B) would not have such surface forces since equation (I196|l 
would then guarantee vanishing toroidal forces everywhere. 

The results found here can be applied directly to precession of neutron stars. For fluid stars, Spitzer (1958) argued that 
precession is inevitable if the magnetic and rotational axes are misaligned; Mestel & Takhar (1972) showed that the star 
precesses about its magnetic symmetry axis with a period Pp = P*/3emagCosx where x is the misalignment angle. For a 
radio pulsar, there would be no effect on the arrival times of pulses if the pulsar beam is along the magnetic axis of the star. 
Wasserman (2003) showed that crustal distortions with a symmetry axis that is also misaligned with the magnetic axis would 
lead to periodically varying timing residuals. For PSR B1828-11, spindown can enhance the effect considerably, and the data 
can be accounted for with B ~ 10^^"^"^ G, x ~ 1 rad, and a modest permanent crustal distortion ~ 0.01 times the magnetic 
distortion. (Perhaps fortuitously, this is close to the crustal distortion found by Cutler et al. 2003 for relaxation near the actual 
rotation frequency of PSR B1828-11.) The model favors prolate figures (see also Akgiin et al. 2006), as would be expected 
from (predominantly) toroidal fields. Why the magnetic and spin axes are misaligned remains unexplained. Moreover, the 
effects of the slow, time variable ffuid motions that would be required in such a model (e.g. Mestel & Takhar 1972; Mestel et 
al. 1981; Nittmann & Wood 1981) have yet to be computed. 

In this paper, we have not examined the effects of rotation, internal velocity fields, multi-fluid components, drag and 
dissipation. These will likely introduce new modes and will alter the properties of modes of non-rotating stars. 
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APPENDIX A: STABILITY CRITERIA FOR A MAGNETIC FIELD H{p, B) 



The coefRcients for m = for a magnetic free energy F{p,B) are given by equation (|108|l . where the various quantities are 
defined in equations (|104p and (|107p . The hydrostatic equihbrium in the absence of magnetic fields is sphericaUy symmetric, 
drPo + podr<l>o = 0. In the presence of magnetic fields, the equilibrium is given by equation (I96|l . Using these equations, we 
can rewrite the coefRcients as, to lowest order in _F, 



bo~T2 



d\np 
dr 
d\np 
dr 



dr 

Ui cot e-U2 + 2U3 cot e 



(Al) 



Co ~ -U2 cot 9 + Us cot e 
where. 

To = B'Fbs + 2pBF^pB + p'F^, 



BF, 



- BdrBF^BB - BdrpF^pB - pdrBF^pB - pdrpF^p 



(A2) 



Ta = 2[/o cot 6 - -{BF,b cot 9 + BdgBF,BB + BdgpF,pB + pdgBF,pB + pdgpF,pp) 



Consider the case of a magnetic field H{p, B). In this case, the magnetic free energy F{p, B) is given through H = AttF^b 
(equation [2]|. To lowest order in F, the density is a function of radius, p{r). Therefore, partial derivatives of p with respect to 
the angle 9 can be dropped. Then, equation (|Aip can be written equivalently as. 



1 

cot 

bo 



cot^9 



din p 
dlnr 

dlnp 
dlnr 



^ din p 
dmr 

Q3 + Q5 



(A3) 



Co 



where, we define. 



BF,B + B^Fj 



)i = B^F^BB + pBF,pB 



Q2 = Qo + Oi 1 



1 - 



dlnB 
d In r 



dlnB 
91nr 



Qi 

1 - 



^ _ dlnB 

d In sin 9 
dlnB 



d In sin 9 



B{r,9) = 



The magnetic induction is given by equation (|15p. which we can rewrite in terms of a new arbitrary function g as, 
pg[Hr sin 9) 



H 



Let C, — Hr sin 9 be the argument of the function g, and define. 



V 



ding 



dlnp 



dlnH 



and 



d in ' d In r ' 9 In p 

Then, after some algebra it follows that, 
dlnB _ ^il-ap)+r]{l+^ap) , dlnB 



CTB — 



dlnH 
dlnB 



dlnr 1 + cr _g(l — ??) 

Using H — 4n F,B, we also get, 

and Qi 



^ HB , 

Qo = — 1 + crs 
471 



and 



KB 

An 



9 In sin l + (7_g(l — ?;) 



(cTp + as) 



(A4) 

(A5) 
(A6) 
(A7) 
(A8) 



Then, the coefficients become (from equation IA3 | 
HB [l + aB+ ^{(Tp + crfl)]2(l - ri) 



fflo 



bo 



CO 



Airr^ 1 + crs (1 — r;) 

HB cot 9[l + aB+ ^{(7p + ^b)]{1 + (tb){1 - v) 
2nr^ 1 + 0-5(1-77) 

HBcot^9 {l + aB f{l - n) 



(A9) 



47^7-2 



1+ (713(1-77) 



Since 60 = 4aoCo, one of the stability conditions is immediately marginally satisfied. The other two conditions give, 
1-77 



1 + 0-3(1-77) 



> 



(AlO) 
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Using equation HA7|) we can rewrite this condition as, for 1 + as > 0, 
dlnB 



d In sin 8 



< 1 . (All) 



Thus, the magnetic fields are marginally stable for B oc sinO. 

For a strongly type II superconducting star H — H[p), so that gb ~ 0, and equation (|A10|I reduces to 77 < 1. For a normal 
conducting star H — B, so that era — 1, and we get (1 — r7)/(2 — r;) > 0. This condition can be expressed in an alternative 
way by noting that equation (|A5|) for a normal conductor is B = pg{Brsm9)/B. Thus, B is given as a function of itself. This 
equation can be rewritten as B = h{pr^ sin^ 6) /r sin 9, and the magnetic free energy is given hy F = B^/8n = p f {pr^ sin^ 0), 
where h and / are arbitrary functions. From here and from equation (|A7|I it follows that, defining zu = pr^ sin^ 0, 

dlnB _ din/ _ ri_ ^^^2) 



d In sin 9 dlnzu 2 — 7; 

The same result is obtained by considering the derivative of B with respect to r, though it involves more algebra. Thus, the 
stability condition for the normal conducting case is better expressed as, 

< 1 . (A13) 

d In ti7 

For a normal conducting star, the field is marginally stable for f oc zu, i.e. B oc pr-sin6', as noted in MA\ Similarly, for a 

strongly type II superconducting star, the field is marginally stable for g oc i.e. B oc prsin9, as noted in M.5\ 



